1W-fll64  834 


UNCLASSIFIED 


SIHPLE  HOOELS  FOR  SURF20NE  SEOIHENT  TRAHSPORT<U>  NAVAL 
CIVIL  ENOINEERINO  LAB  PORT  HUENENE  CA  J  A  BAILARO 
DEC  8S  NCEL-TN-174B 


MICROCOPY  RESOLUTION  TEST  CHART 

•-•’-'tyai  UiiP'*"  Af  <‘t»N0AOn$  A 


NCEL 

Technical  Note 


N-  1740 
December  1985 
By  James  A.  Bailard 

Sponsored  By  Director  of 
Navy  Laboratories 


Simple  Models  for  Surfzone 
Sediment  Transport 


ABSTRACT  Energetics-based  longshore  and  on-offshore  transport  models 
have  been  developed  for  the  surfzone.  The  longshore  transport  model  des¬ 
cribes  the  distribution  of  the  longshore  transport  rate  across  the  surfzone. 
Spatial  integration  of  this  distribution  produces  an  equation  similar  in  form 
to  the  wave  power  equation  but  with  a  wave  power  coefficient  which  is  a 
function  of  the  breaker  angle,  the  beach  slope,  and  the  ratio  of  the  orbital 
velocity  magnitude  at  the  breakpoint  divided  by  the  fall  velocity  of  the 
sediment.  The  on-offshore  transport  model  describes  the  magnitude  and 
direction  of  the  spatially  averaged  on-offshore  sediment  transport  rate  in 
the  surfzone.  These  parameters  were  found  to  be  functions  of  the  incident 
wave  height  and  the  mean  beach  slope.  Both  the  longshore  transport  and 
the  on-offshore  transport  models  have  been  incorporated  into  a  simple 
two-line  shoreline  evolution  model. 


NAVAL  CIVIL  ENGINEERING  LABORATORY  PORT  HUENEME  CALIFORNIA  93043 


METRIC  CONVERSION  FACTORS 


Unclassified 


SECu«i^'f  CLASS<fiC*'^iO»<  O*’  ^MiS  P*Gf  '»>>»«  P«/«  Fnfn^d} 


REPORT  DOCUMENTATION  PAGE 

KFAr>  iNSTRLCnONS 

HFKJRF  C  'lStrj.I-TI.M.  FORM 

1  REPORT  SJM0FB  b  GOVT  ACCESSION  NO 

TN-1740  1  DN387323 

3  RFCiR’EnT'*,  '  a'-Al")'.  n.'MBFR 

*  title  S.ibr.ilr 

SIMPLE  MODELS  FOR  SURFZONE 

SEDIMENT  TRANSPORT 

S  Type  of  report  pfriOO  COvERFO 

Final;  Oct  1983  -  Sep  1984 

6  PERFORMING  ORG  REPORT  NUMBER 

7  AuThQRy 

James  A.  Bailard 

6  contract  o»  grant  NuMBER'ii 

9  PENEORMING  ORGANIZATION  NAME  ANO  AOORE9S 

NAVAL  CIVIL  ENGINEERING  LABORATORY 

Port  Hueneme,  California  93043 

»o  PROGRAM  Element  project  task 

AREA  &  WORK  UNIT  NtiMSCRS 

61152N; 

ZROOO-01-187 

1’  COn’ROcliNG  OPPiCE  NAME  anO  aOORESS 

Director  of  Navy  Laboratories 

Washington,  DC  20360 

b'  REPORT  OATg 

December  1985 

1  j  numbe r  of  Pages 

59 

t«  MCN'ToRinG  agency  name  8  ADDRESS'''/  ditlerent  Irnm  Office) 

15  security  Class  (of  ihi$  report) 

Unclassified 

ISa  declassification  DOWNGRADING 
schedule 

16  OlSTBiBUTiON  STATCmEN’  I-jI  fh«»  Rfporf- 


Approved  for  public  release;  distribution  unlimited. 


1?  OtSTBieijTiO»J  6T*TCmCnT  fo/  fhr  »r>i^rerf  i**  Blorh  ?0.  >/  different  from  fftporfi 


18  S. /BP;.  fMEKTABY  NOTES 


KEY  «OR0S  '  Conrtnu*  oo  rtvrrM*  iid*  if  nee*asmry  »nd  identify  t»v  bferh  number) 

Model,  surfzone,  beach,  sediment  transport,  longshore  transport,  on-offshore  transport, 
line  model 


20  abstract  fConfinue  «n  reeetee  tide  U  neceetery  and  identify  by  Mock  lYtrmberJ 

Energetics-based  longshore  and  on-offshore  transport  models  have  been  developed  for 
the  surfzone.  The  longshore  transport  model  describes  the  distribution  of  the  longshore 
transport  rate  across  the  surfzone.  Spatial  integration  of  this  distribution  produces  an 
equation  similar  in  form  to  the  wave  power  equation  but  with  a  wave  power  coefficient 
which  is  a  function  of  the  breaker  angle,  the  beach  slope,  and  the  ratio  of  the  orbital  velocity 

continued 

00,j?C“,  M73  COITION  or  I  iiiov»j  ij  oaioLtTc  Unclassified 

fCCuUlTr^L A»liriCATtON  OF  THIS  PAG!  0«l«  Entered) 


magnkude  at  the  breakpoint  divided  by  the  fall  velocity  of  the  sediment.  The  on-offshore 
transport  model  describes  the  magnitude  and  direction  of  the  spatially  averaged  on-offshore 
sediment  transport  rate  in  the  surfzone.  These  parameters  were  found  to  be  functions  of  the 
incident  wave  height  and  the  mean  beach  slope.  Both  the  longshore  transport  and  the  on- 
offshore  transport  models  have  been  incorporated  into  a  simple  two-line  shoreline  evolution 
model. 


library  Card 

Naval  Civil  Engineering  Laboratory 

SIMPLE  MODELS  FOR  SURFZONE  SEDIMENT  TRANSPORT 

(Pinal),  by  James  A.  Bailard 

TN’1740  59  pp  illus  December  1985  Unclassified 

I.  Sediment  transport  2.  Line  model  I.  ZR(XX)-01-187 

Energetics-based  longshore  and  on-offshore  transport  models  have  been  developed  for  the 
surfzone.  The  lonj^hore  transport  model  describes  the  distribution  of  the  longshore  transport 
rate  across  the  surfzone.  Spatial  integration  of  this  distribution  produces  an  equation  similar  in 
form  to  the  wave  power  equation  but  with  a  wave  power  coefficient  which  is  a  function  of  the 
breaker  angle,  the  beach  slope,  and  the  ratio  of  the  orbital  velocity  magnitude  at  the  breakpoint 
divided  by  the  fall  velocity  of  the  sediment.  The  on-offshore  transport  model  describes  the  mag¬ 
nitude  and  direction  of  the  spatially  averaged  on-offshore  sediment  transport  rate  in  the  surfzone. 
These  parameters  were  found  to  be  functions  of  the  incident  wave  height  and  the  mean  beach 
slope.  Both  the  longshore  transport  and  the  on-offshore  transport  models  have  been  incorporated 
into  a  simple  two-line  shoreline  evolution  model. 


Unclassified 


SeCuntTy  CCASSiriCATlON  or  this  PAcenrhan  Data  Fntaratf) 


INTRODUCTION 


Quantitative  models  for  nearshore  sediment  transport  processes  are 
necessary  for  understanding  and  predicting  nearshore  bathymetry  changes. 

The  latter  include  both  beach  profile  changes  and  changes  in  the  planform 
shape  of  the  beach.  Changes  in  the  beach  profile  are  most  often  associated 
with  variations  in  the  incident  wave  conditions,  while  planform  changes 
are  associated  with  an  interruption  of  the  littoral  drift  (e.g.,  by  a 
groin  or  harbor  jetty). 

Beach  profile  changes  have  time  scales  which  range  from  a  few  days 
when  associated  with  a  single  storm,  to  several  months  when  associated 
with  seasonal  wave  changes.  Fluctuations  in  the  beach  profile  are 
primarily  the  result  of  on-offshore  sediment  transport  processes  and 
take  place  more  or  less  along  the  entire  coastline. 

Beach  planform  cl.anges,  on  the  other  hand,  have  a  time  scale  which 
varies  with  distance  from  the  littoral  barrier.  For  example,  in  the 
case  of  a  single  groin  the  initial  planform  changes  are  restricted  to 
the  immediate  vicinity  of  the  groin.  Over  time,  the  effects  of  the 
groin  are  felt  at  increasing  distances  downcoast.  Depending  on  the 
overall  size  of  the  groin  and  the  rate  of  longshore  sediment  transport, 
the  initial  effects  are  evident  within  a  few  weeks  while  the  cumulative 
downcoast  effects  sometimes  take  several  years  to  develop. 

The  most  widely  accepted  surfzone  sediment  transport  models  are 
simple,  empirically  based  expressions  which  can  be  readily  applied  to 
many  conunon  coastal  engineering  problems.  Unfortunately,  the  simplicity 
of  these  models  and  their  lack  of  precise  physical  basis  can  lead  to 
erroneous  results  when  they  lack  the  flexibility  and  degree  of  detail 
needed  for  a  number  of  applications.  More  sophisticated  models  have 
been  proposed,  however,  these  models  are  so  complex  that  they  have  had 
only  limited  application.  Moreover,  nearly  all  existing  longshore  and 


on-offshore  sediment  transport  models  lack  a  unity  in  their  development, 
with  the  two  transports  treated  as  distinct  processes  instead  of  different 
components  of  the  same  process. 

The  objective  of  the  present  study  was  to  develop  a  series  of 
analytical  and  numerical  models  to  predict  the  longshore  and  on-offshore 
sediment  transport  rates  as  well  as  the  evolution  of  nearshore  bathymetry. 
The  longshore  and  on-offshore  sediment  transport  model  developments  were 
to  share  a  common  energetics  basis  while  the  bathymetry  evolution  model 
was  to  be  based  on  existing  line  model  methodology.  Although  the  models 
were  to  incorporate  as  much  of  the  physics  of  the  sediment  transport 
processes  as  possible,  the  resulting  equations  were  to  be  simple  enough 
to  be  easily  applied  by  practicing  coastal  engineers.  Funding  for  this 
study  and  earlier  supporting  work  was  provided  by  the  Independent  Research 
Program  of  the  Naval  Civil  Engineering  Laboratory. 

This  report  consists  of  the  following  sections.  First,  existing 
models  for  longshore  transport,  on-offshore  transport,  and  nearshore 
bathymetry  evolution  are  reviewed.  Next,  a  generalized  form  of  Bagnold's 
energetics-based  sediment  transport  model  for  streams  is  reviewed,  and 
then  applied  to  idealized  surfzone  conditions.  The  resulting  equations 
for  the  local  longshore  and  on-offshore  sediment  transport  rat  are 
expressed  as  functions  of  several  local  wave,  current,  and  sediment 
parameters.  These  local  equations  are  simplifed  using  measured  surfzone 
velocity  field  data,  and  combined  with  existing  surfzone  wave  and  current 
models  to  predict  the  spatially  integrated  longshore  and  on-offshore 
sediment  transport  rates  as  a  function  of  incident  wave  parameters. 
Finally,  both  the  longshore  and  the  on-offshore  transport  models  are 
combined  with  existing  line-model  methodology  to  numerically  predict  the 
evolution  of  nearshore  bathymetry  as  a  function  of  time. 


EXISTING  SURFZONE  SEDIMENT  TRANSPORT  MODELS 

Longshore  Transport  Models 

Existing  longshore  transport  models  fall  into  two  general  categories, 
those  that  predict  the  spatially  integrated  longshore  transport  rate, 


and  those  that  predict  the  distribution  of  the  local  longshore  transport 
rate.  The  most  widely  accepted  longshore  transport  model  is  the  wave 
power  equation,  alternatively  called  the  Scripps  equation  or  the  CERC 
equation  (Inman  and  Bagnold,  1963;  U.S.  Army  CERC,  1977).  This  equation 
linearly  relates  the  spatially  integrated  immersed  weight  transport 
rate,  1^^,  to  what  has  been  termed  the  longshore  component  of  wave  energy 
flux,  i  .e.  , 


where  the  wave  power  coefficient,  K,  is  generally  assumed  to  be  a  constant 
equal  to  0.77  (Komar  and  Inman,  1970).  The  factor  can  be  expressed 
as 


P.  =  (E  Cn).  sin  a  cos  a 
Si  bob 

where;  E  =  energy  of  the  incident  waves 
Cn  =  wave  group  velocity 
Of  =  incident  wave  angle 
=  value  at  breaking 


(2) 


The  volumetric  transport  rate, 
weight  transport  rate  by 


Q^,  is  related  to  the  immersed 


=  p  g  (S^-1) 

where:  p  =  fluid  density 

g  =  gravity 

=  specific  gravity  of  the  sediment 


(3) 


relative  to  the  fluid 


N 


at  rest  volume  concentration  of  sediment  on  the  beach  face 


Although  Komar  and  Inman  (1970),  Dean  et  al.  (1983),  Kraus  et.  al. 
(1983)  and  others  have  shown  that  Equation  1  does  an  adequate  job  of 
predicting  the  overall  trend  of  longshore  transport  rate  field  data,  it 
significantly  overpredicts  longshore  transport  rate  for  laboratory  data 
(see  Figure  1).  Moreover,  each  of  the  above  mentioned  field  studies 
suggests  a  somewhat  different  value  for  the  wave  power  coefficient,  K. 

A  possible  explanation  for  the  inability  of  Equation  1  to  predict 
both  laboratory  and  field  measurements  of  longshore  transport  rates 
using  a  single  value  of  K  is  that  sediment  transport  in  the  laboratory 
may  he  mostly  bedload  while  transport  in  the  field  may  be  mostly  suspended 
1  oad  . 

The  question  of  the  relative  importance  of  bedload  versus  suspended 
load  transport  is  a  subject  of  open  debate.  A  number  of  investigators 
(e.g.,  Komar,  1976)  assert  that  most  of  the  transport  in  the  field  is  by 
bedload.  In  truth,  this  question  may  be  one  of  semantics,  for  while 
bedload  is  precisely  defined  as  all  sediment  which  is  supported  in 
transport  by  grain  to  grain  collisions,  in  practice  bedload  has  often 
been  defined  as  all  sediment  transported  within  10  cm  of  the  bottom. 
Obviously  this  latter  definition  includes  sediments  transported  both  as 
bedload  and  as  suspended  load. 

The  wave  power  equation  predicts  the  spatially  integrated  longshore 
transport  rate.  A  number  of  longshore  transport  models  have  been  developed 
to  predict  the  distribution  of  the  longshore  sediment  transport  across 
the  surfzone.  Nearly  all  of  these  models  utilize  an  adaptation  of  a 
stream-based  sediment  transport  model.  For  example,  the  models  by 
Ostendorf  and  Madsen  (197)9  and  Swart  (1976)  utilize  traction-based 
models,  while  the  models  by  Komar  (1977),  Thornton  (1973),  and  Bailard 
(1981)  utilize  energetics-based  models.  In  order  to  predict  the 
distribution  of  the  longshore  transport  rate  across  the  surfzone,  these 
models  incorporate  various  methods  for  predicting  wave-induced  velocities 
and  the  longshore  current.  Because  of  this  added  complexity,  most  of 
the  models  require  numerical  solution  and  have  not  enjoyed  wide  spread 
application. 


4 


On-offshore  Sediment  Transport  Models 

Unlike  longshore  transport  models,  quantitative  models  for  predicting 
the  on-offshore  sediment  transport  rate  on  a  beach  are  scarce.  One 
difficulty  is  that  the  time-averaged  on-offshore  sediment  transport  rate 
represents  a  small  difference  between  the  two  large  instantaneous  onshore 
and  offshore  transports  induced  by  each  wave.  As  a  result,  secondary 
features  in  the  flow  field,  such  as  wave  velocity  asymmetries  and  weak 
on-offshore  currents,  have  a  far  greater  effect  on  the  on-offshore 
sediment  transport  rate  than  on  the  longshore  sediment  transport  rate. 

Conceptually,  the  on-offshore  sediment  transport  in  the  surfzone 
can  be  viewed  as  a  balance  between  the  downslope  component  of  gravity, 
the  asymmetry  of  the  orbital  velocity  field,  and  the  time-averaged 
offshore  current.  Although  the  above  factors  (and  others)  have  been 
recognized  for  some  time,  the  complexity  of  the  fluid/sediment  motions 
in  the  surfzone  have  forced  many  investigators  to  seek  more  qualitative 
relationships  between  the  incident  wave  properties  and  the  direction  of 
on-offshore  sediment  transport.  Examples  of  this  approach  include  the 
wave  height  models  by  Saville  (1957),  Aubrey  (1978),  and  Short  (1978), 
and  the  wave  steepness  models  by  Dean  (1973)  and  Hattori  and  Kawamoto 
(1981). 

A  few  quantitative  models  have  been  proposed  including  those  by 
Swart  (1974),  Bowen  (1980)  and  Bailard  (1981).  The  latter  two  models 
are  based  on  the  energetics  concepts  of  Bagnold  (1963,  1966)  which  are 
discussed  in  the  next  section.  Swart's  on-offshore  transport  model  is 
based  on  an  extension  of  an  equation  proposed  by  Bakker  (1968)  and 
Bakker  et  al.  (1970).  This  equation  relates  the  on-offshore  sediment 
transport  rate  per  unit  width  of  beach,  q^,  to  the  distance  between  two 
depth  contours,  x^  and  x^ ,  as  follows 

=  Coff  (x^-  x^  +  w)  (4) 

where:  Coff  =  rate  coefficient  assumed  constant  for  a  specific  set  of 

boundary  conditions 


w  =  equilibrium  value  of  (x  “  x  ) 


In  Bakker's  original  formulation  (used  in  a  line  model  development), 
Coff  and  w  were  assumed  to  be  constants,  however,  in  Swart's  model  both 
of  these  parameters  are  allowed  to  vary  as  complex  empirical  functions 
of  the  incident  wave  conditions,  the  local  depth,  the  local  orbital 
velocity  magnitude,  and  the  sediment  size.  The  resulting  model  can  be 
used  to  predict  the  change  in  the  beach  profile  as  a  function  of  time, 
as  well  as  the  shape  of  the  equilibrium  beach  profile.  The  complexity 
of  Swart's  model,  however,  has  prevented  its  widespread  application. 

Line  Models 

Line  models,  which  represent  the  nearshore  bathymetry  by  one  or 
more  depth  contour  lines,  are  a  simple  method  of  predicting  nearshore 
bathymetry  changes  as  a  function  of  time.  In  a  one-line  model,  the 
entire  beach  is  represented  by  the  shoreline  (MSL)  contour  (see  Figure  2). 
The  beach  profile  is  assumed  constant,  on-offshore  sand  movements  are 
neglected,  and  all  shoreline  changes  are  associated  with  the  divergence 
of  the  longshore  transport  rate. 

Pelnard  Considere  (1956)  developed  the  first  one-line  model  by 
combining  a  linearized  longshore  transport  equation  with  the  continuity 
equation  for  sand.  The  resulting  equation  for  the  shoreline  position, 

X,  took  the  form  of  the  diffusion  equation. 

a  9^2 

^  -  o  9  X  ,  . 

at  D  2  ^  ^ 

9y 

where:  x  =  longshore  coordinate 

q  =  derivative  of  the  longshore  transport  rate  with  respect 
to  the  angle  of  wave  incidence 

D  =  depth  out  to  which  sand  transport  is  negligible 

Equation  5  has  been  solved  for  a  variety  of  boundary  and  initial 
conditions  as  detailed  in  Bakker  (1968),  Bakker  et  al.  (1970),  and 
LeMehaute  and  Soldate  (1977). 


Multiple  line  models,  depicting  two  or  more  depth  contours,  provide 
information  on  the  local  beach  slope  (or  profile)  and  provide  a  mechanism 
for  incorporating  on-offshore  sediment  movements.  Bakker  (1968)  developed 
the  first  two-line  model  where  the  coupled  equations  governing  the 
position  of  each  line  were  as  follows 


(V°i)  aT 


°2  at  ay  ‘^x 

where  (referring  to  Figure  3) 


Xj  =  distance  of  the  first  (MSL)  contour  line  from  the 
base  line 

X2  =  distance  of  the  second  (D^)  contour  line  from  the  base 
line 

Dj  =  depth  associated  with  the  second  contour  line 


)g  =  height  of  the  sand  berm  above 


Dj  +  =  depth  below  which  sand  transport  is  negligible 

Qj  =  longshore  sediment  transport  rate  associated  with 
the  Xj  contour  line 

=  longshore  sediment  transport  rate  associated  with  the 
contour  line 

=  offshore  sediment  transport  rate  from  the  Xj  to  the 
^  contour  line 


Solution  of  Bakker 's  two-line  model  is  most  easily  accomplished 
numerically.  The  resulting  solutions  correctly  simulate  local  beach 
bathymetry  changes  which  occur  around  a  groin,  including  the  steepening 
of  the  profile  on  the  upcoast  side  and  the  flattening  of  the  profile  on 
the  downcoast  side  (see  Figure  4). 


Recently,  Perlin  and  Dean  (1983)  have  developed  a  six-line  model 
which  is  a  logical  extension  of  Bakker’s  two-line  model.  Although  the 
six-line  model  more  accurately  represents  the  shape  of  the  beach  profile 
than  the  two-line  model,  neither  model  accurately  predicts  the  beach 
profile  changes  associated  with  changes  in  the  incident  wave  field. 

ENERGETICS  SURFZONE  SEDIMENT  TRANSPORT  EQUATIONS 

Because  of  the  complexity  of  the  flow  field  in  the  surfzone  and  the 
simplicity  of  Bagnold's  energetics-based  sediment  transport  concepts, 
the  latter  has  found  widespread  application  in  modeling  surfzone  sediment 
transport  (Inman  and  Bagnold,  1963;  Komar,  1977;  Thornton,  1973;  Bowen, 
1980;  Bailard,  1981).  The  following  is  a  brief  summary  of  Bagnold's 
(1963,  1966)  sediment  transport  model  for  streams,  and  Bailard's  (1981, 
1983)  application  of  this  model  to  sediment  transport  in  the  surfzone. 

Generalized  Transport  Equation 

Bagnold's  energetics-based  sediment  transport  model  for  streams 
assumes  that  the  transport  occurs  as  two  distinct  modes,  bedload  and 
suspended  load.  Bedload  sediment  transport  occurs  as  a  thin  granular- 
fluid  layer  which  is  supported  by  the  bed  via  grain  to  grain  collisions. 
Suspended  load  transport  takes  place  above  the  bedload  layer,  where  the 
sediment  is  supported  by  fluid  turbulence.  In  both  cases,  energy  is 
expended  by  the  stream  in  transporting  sediment,  thus  an  analogy  can  be 
made  of  a  stream  acting  as  a  sediment  transport  machine.  Bagnold  assumed 
that  the  rate  of  sediment  transport  was  proportional  to  the  rate  of 
energy  dissipation  of  the  stream.  The  proportionality  constants  associated 
with  the  bedload  and  suspended  load  transports  were  termed  efficiency 
factors,  assumed  to  be  constant.  Bailard  (1981)  generalized  Bagnold's 
energetics-based  stream  model  for  time-varying  flow  over  an  arbitrarily 
sloping  bottom.  The  resulting  sediment  transport  equation  was 
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where:  =  instantaneous  sediment  transport  rate  vector 

p  =  density  of  water 

=  drag  coefficient  of  the  bed 
£  =  bedload  efficiency  factor 

D 

<|i  =  internal  angle  of  friction  for  the  sediment 
tan  p  =  bed  slope 

=  suspended  load  efficiency  factor 
W  =  fall  velocity  of  the  sediment 

=  instantaneous  nearboi.tom  fluid  velocity  vector 
i  =  unit  vector  directed  upslope 
<  >  =  a  time-average  quantity 


An  important  feature  of  the  above  equation  is  that  the  bedload 
transport  (first  bracketed  quantity)  and  the  suspended  load  transport 
(second  bracketed  quantity)  consist  of  primary  components  directed 
parallel  to  the  instantaneous  fluid  velocity  vector,  and  secondary 
components  directed  downslope.  The  latter  are  associated  with  the 
downs  lope  component  of  the  sediment  load. 


Local  Surfzone  Transport  Equations 

Figure  5  depicts  a  plane  contour  beach  with  the  x-axis  directed 
shoreward  and  normal  to  the  beach,  and  the  y  axis  directed  parallel  to 
the  beach.  The  slope  of  the  beach  is  tan  p.  An  incident  wave  with  deep 
water  angle,  a  ,  is  assumed  to  have  a  local  wave  angle  a. 


For  simplicity,  the  nearbottom  velocity  field  is  assumed  to  be 
composed  of  an  oscillatory  component  u  oriented  at  an  angle  a  to  the 
x-axis,  and  steady  velocity  components  u  and  v  directed  in  the  x  and  y 
directions,  respectively.  The  total  velocity  vector  becomes 


(u  cos  a  +  u)  i  +  (u  sin  a  +  v)  j 


(9) 


In  addition,  the  oscillatory  velocity  component  is  assumed  to  be 
composed  of  a  primary  component  u^  with  frequency  a  and  higher  harmonics 
u^2  with  frequencies  2a,  etc. ,  so  that 

u  =  u  cos  ot  +  u  „  cos  2  at  +  . . .  (10) 

m  mz 


Substitution  of  Equation  9  into  Equation  8,  and  assuming  that 
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then  the  idealized  on-offshore  and  longshore  transport  equations  become 
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Equations  11  and  12  describe  the  local  on-offshore  and  longshore 
transport  rates  as  a  function  of  local  wave,  current,  beach,  and  sediment 
parameters.  Combined  with  nearbottom  current  meter  records,  these 
equations  can  be  used  to  predict  the  local  sediment  transport  rate 
vector  as  a  function  of  time.  For  modeling  purposes,  however,  the 
distribution  and  spatially  integrated  values  of  the  time-averaged 
sediment  transport  rate  are  of  greater  interest.  Predicting  these 
transport  rates  as  a  function  of  the  incident  wave  field  requires  the 
use  of  additional  surfzone  model  theory. 

Considering  the  longshore  transport  Equation  12,  all  of  the  free 
parameters  can  be  predicted  from  linear  wave  shoaling/ ref raction  theory 
and  one  of  several  radiation  stress-based  longshore  current  models  which 
are  available  (e.g.,  Longuett-Higgins ,  1970).  In  its  full  complexity, 
this  procedure  results  in  a  complex  expression  for  the  longshore  sediment 
transport  rate  which  must  be  numerically  evaluated  (Bailard,  1981).  As 
a  result,  the  model  is  difficult  to  apply  and  the  results  are  not  easily 
generalized.  In  the  next  section,  the  approach  used  by  Bailard  (1981) 
will  be  simplified  yielding  a  less  complex  and  more  easily  applied 
longshore  transport  model. 


Equation  11,  which  predicts  the  on-offshore  transport,  is  more 
difficult  to  apply  than  the  longshore  transport  Equation  12  because  the 
wave  velocity  moments,  4*^  and  *^he  mean  velocity,  6^,  cannot  be 

predicted  from  existing  surfzone  model  theory.  The  approach  taken  in 
the  present  study  was  to  estimate  these  parameters  from  measured  field 
data  and  correlate  their  values  with  the  incident  wave  characteristics. 
Details  of  this  procedure  will  be  discussed  in  a  later  section. 


SPATIALLY  INTEGRATED  LONGSHORE  TRANSPORT 

The  ubiquity  of  the  wave  power  equation  (Equation  1)  suggested  that 
the  present  longshore  transport  model  should  take  the  form  of  a  modifi¬ 
cation  to  the  wave  power  equation.  Following  the  approach  by  Bailard 
(1981),  an  expression  was  sought  for  the  wave  power  coefficient,  K,  as  a 
function  of  incident  wave,  beach,  and  sediment  parameters. 

The  principal  obstacles  to  obtaining  a  close-form  solution  for  K 

*  * 

were  the  complex  expressions  for  the  wave  velocity  moments,  and  u^. 

Guza  and  Thornton  (1981)  have  shown,  however,  that  for  weak  mean  currents 

*  * 

and  a  Gaussian  orbital  velocity  distribution,  u  and  u  are  constants 

^  ^  *  * 

equal  to  0.58  and  1.13,  respectively.  Field  measurements  of  u^  and  u^ 
(Bailard,  1983)  support  this  finding  with  measured  values  equal  to  0.60 
to  1.25,  respectively.  Consequently,  an  important  simplifying  assumption 
made  in  the  present  study  was  that  u^  and  u^  are  constants  equal  to  the 
above  measured  values. 

The  longshore  current  model  used  in  the  present  development  was  the 
radiation  stress-based  model  by  Ostendorf  and  Madsen  (1979).  In  order 
to  obtain  a  tractable  solution,  a  number  of  simplifications  were  made 
with  regard  to  the  model's  provisions  for  finite  breaker  heights  and 
finite  longshore  current  strengths.  The  resulting  expression  for  the 
longshore  current,  v,  was  as  follows 


V  =  V  V 

c 


« 
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where  is  a  characteristic  longshore  current  strength  and  v-  is  the 
dimensionless  longshore  current  distribution  expressed  as  follows  (P  ^  0.4) 


/I  -  C  \  3 


The  parameters  C^,  C^,  and  are  functions  of  the  lateral  mixing 


parameter,  P,  as  follows 
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and  the  dimensionless  surfzone  position,  x^,  is  defined  as 


X,  -  X 

b  s 


where  x^  and  x^  are  the  positions  of  the  shoreline  and  the  break-point 
respectively . 

The  characteristic  longshore  current  strength,  v^,  is  defined  as 


V  =  6  u  . 

c  c  mb 


where  u^^  is  the  oscillatory  velocity  magnitude  at  the  breakpoint 
obtained  from  shallow  water  wave  theory,  i.e., 


and  6  is  the 
c 


relative  longshore  current  strength  parameter  defined  as 
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5n  tan  A  sin  a. 
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(19) 


Note  that  in  Equations  18  and  19,  tan  A  is  equal  to  the  average  beach 
slope,  tan  p,  modified  for  wave  setup,  i.e.. 


tan  A 


tan  p 
,  ^  3  2 

^  8 


(20) 


and  is  the  ratio  of  the  surfzone  wave  height,  H,  to  water  depth,  h, 
defined  as 


y  =  Yf,  X  <1 

>  =  \{—)  * 

and  is  assumed  to  be  equal  to  0,8. 

Combining  the  longshore  current  Equation  13  and  the  local  longshore 
transport  Equation  12,  and  assuming  that  Snells  Law  applies,  the  following 
equation  is  obtained  for  the  cross-shore  distribution  of  the  longshore 
transport  rate 
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The  spatially  integrated  longshore  transport  rate, 

by  integrating  Equation  19, 

may  be  obtained 
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Noting  that  K  =  tan  =  sin  and 
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P  tanA  ‘^mb  8  ^b  ^ SL 

(24) 

then  the  wave  power  coefficient  K  can  be  expressed  as 

K  =  e„  K  +  K„  +  K, 

(25) 
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Further  simplifications  to  Equations  26  through  28  are  possible 
using  the  results  of  Komar  (1976).  Komar  analyzed  the  results  of  labo¬ 
ratory  and  field  measurements  of  longshore  currents  and  found  that  P  was 
approximately  equal  to  0.2  and  the  ratio  tan  p/c^  was  approximately 
equal  to  7.  Introducing  these  assumptions  into  Equations  26  through  28 
the  following  results  are  obtained 


Kj  =  0.385  +  20  sin  20^^  +  0.074  tan  p 


=  0.228  -55- 


K^  =  0.123  tan  P 


(¥) 
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The  only  free  parameters  remaining  in  Equation  25  that  must  be 

specified  are  the  bedload  and  suspended  load  efficiency  factors.  Bagnold 

(1966)  found  that  Eg  and  Eg  were  equal  to  0.15  and  0.01  respectively  for 

stream  flow.  In  the  present  case,  these  factors  were  estimated  from 

laboratory  and  field  measurements  of  the  longshore  transport  rate. 

The  laboratory  data  used  to  estimate  e„  and  e  were  selected  data 

D  S 

from  Saville  (1949,  1950)  and  from  Shay  and  Johnson  (1951).  The  field 
data  were  from  Bruno  et  al.  (1980),  Dean  et  al.  (1983),  Komar  and  Inman 


p 


5.-- 
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(1970),  Kraus  et  al.  (1983),  and  Moore  and  Cole  (1960).  Average 
characteristics  of  these  data  sets  are  summarized  in  Table  1 . 


Initially,  estimates  for  £„  and  e  were  sought  using  a  nonlinear 

D  o 

least  squares  estimation  procedure.  It  was  quickly  found  that  the 
contribution  from  was  negligible.  Eliminating  from  Equation  22 
yields  an  equation  linear  in  £_  and  £  ,  thus  a  simpler  multiple  linear 

D  O 

regression  analysis  was  used.  The  resulting  estimated  values  for  £g  and 
£g,  along  with  their  95%  confidence  intervals,  were 

£„  =  0. 13  ±  0.009 

B 

£g  =  0.032  ±  0.004  (32) 

Utilizing  these  values  in  Equation  25  the  final  equation  for  the  wave 
power  coefficient,  K,  becomes 

9  “mb 

K  =  0.050  +  2.6  sin  2a^  +  .0096  tan  p  +  0.0073  (33) 

- - - - - ^ 

Bedload  Suspended  load 

Figure  6  shows  a  plot  of  the  observed  values  of  K  versus  the 
estimated  values  of  K.  With  the  exception  of  the  El  Moreno  Beach  data 
from  Komar  and  Inman  (1970)  and  the  single  Moore  and  Cole  (1960)  data 
point,  the  fit  is  relatively  good. 

Spatially  Averaged  On-offshore  Transport  Equation 

Equation  11  predicts  the  local  on-offshore  sediment  transport  rate 
as  a  function  of  the  local  wave,  current,  beach,  and  sediment  parameters 
The  steady  onshore  current  is  characterized  by  6^,  while  the  waves  are 
characterized  by  the  orbital  velocity  magnitude,  u  ,  the  total  velocity 

~Jr 

moments,  u^  and  u^ ,  and  the  two  skewness  parameters,  'Hj  and  The 

latter  describe  the  asymmetry  of  the  orbital  velocity  field  as  characterized 


by  the  brief  but  intense  onshore  flow  occurring  beneath  the  wave  crest 
and  the  longer  but  less  intense  return  flow  occurring  beneath  the  wave 
trough . 

In  principal,  u„ ,  u,.  ,  6  ,  and  u  can  be  predicted  from  existing 
3  5  u  m 

surfzone  wave  theory,  while  and  cannot.  For  example,  u^  and  u^ 

have  been  shown  to  be  approximately  constant,  while  U  may  be  assumed  to 

m 

be  a  function  of  the  local  water  depth  (Equation  18).  Moreover,  several 
models  have  recently  been  proposed  to  describe  6^  as  a  function  of  the 
incident  wave  parameters  (e.g.  Svendsen,  1984).  In  practice,  application 
of  these  models  to  field  conditions  is  difficult.  For  example,  field 
measurements  (Guza  and  Thornton,  1981)  show  u  to  be  nearly  constant 

i72 

across  the  surfzone  instead  of  decaying  as  h  as  suggested  by  the 
breaking  wave  theory.  In  addition,  the  models  for  6^  are  too  complex 
for  easy  application  within  the  context  of  developing  a  simple  on-offshore 
transport  model.  As  a  result,  the  approach  taken  in  the  present  stud’'^ 
was  to  empirically  relate  measured  values  of  u  ,  6  ,  f,  ,  and  I*,  wrth 
incident  wave  properties. 

The  field  data  sets  used  to  estimate  the  above  parameters  were 
obtained  from  the  Nearshore  Sediment  Transport  Study  (NSTS)  experiments 
at  Torrey  Pines  Beach,  Calif.,  and  Leadbetter  Beach,  Calif.  (Gable  1979, 
1980).  Each  experiment  consisted  in  part  of  the  simultaneous  measurements 
of  incident  wave  characteristics,  near  bottom  surfzone  velocities,  and 
beach  profile  changes.  The  surfzone  velocities  were  measured  using 
electromagnetic  current  meters  arranged  in  a  cross-shaped  pattern.  One 
line  of  current  meters  was  aligned  parallel  to  the  beach,  while  another 
line  of  current  meters  was  aligned  perpendicular  to  the  beach.  In  the 
present  study,  only  data  from  the  shore-perpendicular  current  meters 
were  used. 

The  data  sets  from  the  two  beaches  differed  significantly  in  terms 
of  the  beach  slope  (0.02  for  Torrey  Pines  Beach  and  0.038  for  Leadbetter 
Beach)  and  the  incident  wave  characteristics.  Figure  7  shows  that  at 
Torrey  Pines  Beach,  the  significant  wave  height  varied  between  0.8  and 
2.0  meters,  while  the  significant  wave  period  varied  between  9  and 
12  seconds.  At  Leadbetter  Beach,  the  significant  wave  height  varied 
between  0.2  and  0.8  meter,  while  the  significant  wave  period  varied 
between  5  and  20  seconds. 


Initially,  it  was  hoped  that  a  set  of  empirical  equations  could  be 

developed  to  predict  u  ,  6  ,  *1*  .  and  as  a  function  of  surfzone  position, 

m  u  1  2  e  . 

incident  wave  properties,  and  the  average  beach  slope.  In  this  way,  the 
distribution  of  the  on-offshore  sediment  transport  could  be  predicted  as 
a  function  of  the  incident  wave  characteristics  for  a  particular  beach 
slope . 

In  order  to  determine  whether  or  not  a  similarity  type  distribution 

existed,  data  was  plotted  as  a  function  of  h/L  ,  x  ,  and  the  Ursell 
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number,  H/k^h  ,  where  H  is  the  wave  height,  h  is  the  water  depth,  is 

the  deep  water  wave  length,  and  is  the  deep  water  wave  number.  None 

of  the  plots  showed  a  significantly  reduced  degree  of  scatter  in  the 

data  relative  to  plots  of  these  parameters  versus  the  depth,  h. 

Because  of  the  poor  results  of  the  above  procedure,  a  less  ambitious 

approach  was  taken  whereby  the  spatially  averaged  values  of  u^,  6^,  4'^ , 

and  4*2  were  correlated  with  gross  properties  of  the  incident  waves.  The 

latter  included  the  incident  significant  wave  height,  incident  significant 

wave  steepness,  and  the  surfzone  similarity  parameter  defined  as 


tan  P 


(34) 


Although  the  results  were  not  conclusive,  the  best  correlation  was 
obtained  when  the  variables  were  plotted  against  the  incident  significant 
wave  height. 

Figures  8  through  11  show  plots  of  u  ,  6  ,  ,  and  versus 
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significant  wave  height  for  each  data  set.  These  plots  suggested  that 
each  of  the  surfzone  parameters  could  be  fitted  to  a  regression  equation 
of  the  following  form 

parameter  =  (a  tan  p  +  b)  +  (c  tan  P  +  d)  (35) 

where:  =  signficant  wave  height,  meters 


tan  p  =  average  beach  slope  within  the  surfzone 


Expanding  F.quation  35  and  using  a  multiple  linear  regression  analysis, 
estimates  for  the  coefficients  a  through  d  were  obtained  for  each  variable 
The  results  are  listed  in  Table  2,  while  the  linear  regression  curves 
are  shown  in  Figures  7  through  10. 

The  regression  equations  for  the  parameters  u  ,  6  ,  .  and  'E  can 

m  u  1  2 

be  combined  with  the  on-offshore  sediment  transport  equation  (Equation  11) 
to  predict  the  spatially-averaged  on-offshore  sediment  transport  rate  in 
the  surfzone  as  a  function  of  the  incident  wave  height  and  the  average 
beach  slope. 

Figure  12  shows  a  plot  of  the  onshore  transport  rate  as  a  function 
ot  the  incident  significant  wave  height  for  three  different  beach  slopes. 
Free  parameter  inputs  used  to  generate  this  figure  were 

=  0.005 

W  =  0.04  m/sec 

S  =2.65 
s 

N  =0.63 
o 

Sb  =  0.13 
=  0.032 

The  behavior  of  the  on-offshore  sediment  transport  model  as  depicted 
in  Figure  11  shows  several  realistic  features.  First,  gradual  onshore 
sand  movement  occurs  when  the  waves  are  small,  while  rapid  offshore 
movement  occurs  with  larger  waves.  This  feature  helps  to  explain  why 
beaches  can  erode  rapidly  under  the  attack  of  large  winter  waves.  In 
contrast,  beaches  accrete  slowly  under  the  action  of  small  summer  waves. 

A  second  feature  of  the  model  is  that  an  equilibrium  wave  height 
exists  which  is  a  function  of  the  average  beach  slope.  With  steep 
beaches,  the  equilibrium  wave  height  is  small,  while  with  flat  beaches 
the  equilibrium  wave  height  is  larger.  A  corollary  to  this  finding  is 
that  given  sufficient  time,  a  beach  will  seek  a  beach  slope  in  equilibrium 
with  the  incident  wave  height.  Figure  13  shows  a  plot  of  the  equilibrium 
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beach  slope  as  a  function  of  the  incident  significant  wave  height. 

Small  waves  produce  a  steep  beach,  while  large  waves  produce  a  flat 
beach.  This  is  in  keeping  with  field  observations. 

One  feature  which  is  lacking  in  the  model  is  any  significant 
variation  in  the  equilibrium  beach  slope  with  sediment  fall  velocity 
(i.e.,  sediment  size).  Varying  the  sediment  fall  velocity  by  an  order 
of  magnitude  has  the  effect  of  changing  the  magnitude  of  the  on-offshore 
transport  rate  (i.e.,  the  rate  at  which  the  beach  responds  to  a  change 
in  wave  height),  but  not  the  equilibrium  wave  height.  Long-term  field 
observations  (e.g.,  U . S .  Army  CERC,  1977)  have  shown  a  link  between 
sediment  size  and  beach  slope,  however,  this  may  be  due  to  threshold  of 
motion  or  percolation  effects,  neither  of  which  are  included  in  the 
present  model . 

It  is  of  interest  to  compare  the  Bakker/Swart  on-offshore  sediment 
transport  equation  (Equation  4)  with  the  present  model.  One  way  to 
accomplish  this  is  to  assume  normally  incident  waves  and  utilize  a 
degenerative  form  of  the  two-line  model  equations.  Assuming  and 
are  zero,  Equations  5  and  6  become 
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Solutions  to  the  above  equations  are  expressed  in  terms  of  the  position 
of  the  contour  lines  and  x^ ,  as  a  function  of  time.  For  present 
purposes,  it  is  more  convenient  to  express  the  solution  in  terms  of  the 
time  history  of  the  average  beach  slope  expressed  as 


Utilizing  the  present  model,  Equations  36  and  37  can  be  combined  with 
the  on-offshore  transport  model  described  by  Equations  11  and  35. 
Assuming  values  for  the  initial  beach  slope,  the  sediment  fall  velocity, 
the  beach  depths  Dj,  D^,  D^,  and  the  drag  coefficient  for  the  bed,  C^, 
these  equations  can  be  time-stepped  ahead  to  predict  the  evolution  of 
the  beach  slope  subject  to  specific  incident  wave  conditions. 

Figures  14  and  15  show  the  results  of  this  procedure  utilizing 
three  different  initial  beach  slopes  and  two  different  incident  wave 
heights.  Input  conditions  for  these  figures  were  as  follows 
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Figure  14 


Figure  15 


0.03;  0.02;  0.015  0.03;  0.02;  0.010 


0.9  meter 
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See  Equation  (45) 
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Note  that  the  equilibrium  beach  slopes  corresponding  to  the  wave  heights 
0.9  and  0.5  meter  are  0.02  and  0.03  respectively.  Generalizing  these 
results,  it  is  seen  that  the  flattening  of  an  oversteep  beach  occurs 
more  quickly  than  the  steepening  of  an  overflattened  beach.  The  rate  of 
change  of  the  beach  slope  at  a  particular  moment  depends  on  both  the 
initial  beach  slope  and  the  incident  wave  height. 

Solution  of  Equations  36  and  38  using  the  Bakker/Swart  on-offshore 
transport  model  (Equation  4)  produces  an  analytic  expression  for  the 
evolution  of  the  beach  slope  as  a  function  of  time  (see  Appendix  A). 

The  resulting  expression  for  the  beach  slope  is 


tan  B.  tan  8 


+  tan  p  e 


where : 


Coff  Coff 


tan  =  initial  beach  slope 

tan  =  equilibrium  beach  slope 

Using  the  same  input  conditions  discussed  above,  except  that  the 
equilibrium  beach  slope  is  specified  as  opposed  to  the  incident  wave 
height  and  Coff  =  0.263  m/day,  the  resulting  time  history  of  the  beach 
slope  is  plotted  in  Figures  14  and  15. 

Comparing  these  results  with  the  results  presented  in  Figures  13 
and  14,  several  differences  are  evident.  First,  in  the  Bakker/Swart 
model,  the  beach  slope  response  curve  is  symmetrical  with  regards  to 
erosion  or  accretion.  The  rate  at  which  the  beach  slope  changes  is  a 
direct  function  of  the  parameter  Coff,  and  the  depths  ,  D^,  and 
Bakker  (1970)  suggested  that  Coff  should  be  between  0.00274  and  0.0274 
m/day,  while  Berlin  and  Dean  (1983)  suggested  a  value  closer  to 
0.263  m/day. 

The  above  results  suggest  that  the  present  model  conforms  more 
closely  to  field  observations  of  beach  profile  response  than  does 
Bakker' s  (1970)  model.  The  latter  could  be  made  more  realistic  by 
allowing  Coff  to  be  a  function  of  the  incident  wave  height,  the  fall 
velocity  of  the  sediment,  and  the  initial  beach  slope.  This  procedure 
was  in  fact  used  by  Swart  (1974),  except  that  his  model  describes  the 
evolution  of  the  complete  beach  profile. 


LINE  MODEL  APPLICATION 


The  present  longshore  and  on-offshore  transport  models  were 
incorporated  into  a  simple  two-line  model  to  predict  the  evolution  of 
nearshore  bathymetry  around  a  single  groin.  Figure  3  shows  the  important 
elements  of  the  line  model.  The  shoreline  (or  mean  sea  level)  contour 
is  described  by  line,  while  the  depth  contour  is  described  by  the 


line.  The  depth,  is  assumed  to  be  equal  to  the  depth  at 

which  active  sand  transport  ceases.  The  elevation  D_  describes  the 

D 

height  of  the  beach  berm  line  above  mean  sea  level. 

The  basic  relationships  governing  the  line  model  include  the 
sediment  continuity  equation,  the  longshore  transport  rate,  and  the 
on-offshore  transport  rate.  Additional  details  discussed  below  include 
breaking  wave  assumptions,  the  depth  of  active  sand  transport,  boundary 
condition,  and  the  development  of  the  controlling  finite  difference 
equations . 


Basic  Equations 


The  sand  continuity  equations  (Equations  5  and  6)  describe  the 
change  in  shore  line  position  as  a  function  of  the  on-offshore  transport 
rate  and  the  divergence  of  the  longshore  transport  rate.  In  these 
equations,  the  on-offshore  transport  represents  an  exchange  of  sand 
across  the  Dj  contour  line. 

The  on-offshore  transport  model  developed  in  the  previous  section 
describes  the  spatially-averaged  on-offshore  sand  transport  in  the 
surfzone.  No  details  regarding  the  distribution  of  this  transport  are 
predicted,  nor  is  a  precise  definition  given  for  the  region  of  application 
in  the  surfzone.  Thus,  for  simplicity  it  was  assumed  that  the  predicted 
on-offshore  transport  rate  was  equal  to  that  crossing  the  depth 
contour  (i.e.,  just  that  required  by  the  two-line  model). 

In  a  two-line  model,  the  longshore  transport  rate  is  divided  into 
two  parts,  and  Q^.  Qj  includes  all  of  the  longshore  transport 
occurring  shoreward  of  the  Dj  depth  contour,  while  includes  all  the 
transport  occurring  seaward  of  this  contour.  As  a  result,  this  division 
of  the  longshore  transport  rate  can  be  expressed  as  follows 


"  "^l  ^2  ”  (I'll)  Qjj 

where:  Q-  =  total  longshore  transport  rate  given  by  Equations  1, 
^  3,  and  33 


(41) 


The  parameter  r|j  is  defined  as 
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where:  x,  =  Dl/h. 
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An  analytical  expression  for  rjj  i®  derived  in  Appendix  B. 


Breaking  Wave  Assumptions 

Wave  input  conditions  for  the  line  model  consist  of  the  incident 
significant  wave  height,  the  significant  wave  period,  and  the  deep  water 
incident  wave  angle  (measured  counterclockwise  from  the  y  axis).  As  the 
waves  progress  into  shallow  water,  refraction  and  shoaling  processes 
cause  the  wave  angle  and  the  wave  height  to  change.  Wave  diffraction  is 
also  important  in  the  vicinity  of  the  groin.  LeMehaute  and  Soldate 
(1980)  and  Kraus  and  Harikai  (1983)  have  developed  methods  for  incorpo¬ 
rating  wave  diffraction  effects  into  single-line  model  methodology.  In 
the  present  model  development,  however,  wave  diffraction  effects  were 
judged  to  be  a  complication  which  could  be  added  later. 

As  the  waves  propagate  into  shallower  and  shallower  water,  the 
ratio  of  wave  height  to  water  depth  eventually  becomes  so  large  that  the 
waves  break.  Following  the  development  of  LeMehaute  and  Soldate  (1980), 
the  breaker  height  and  breaker  angle  are  computed  from  the  deep  water 
wave  height,  H^,  and  the  deep  water  wave  angle  (relative  to  the  x  axis), 
a,  in  the  following  manner 


»b  = 


=  K  K  H 
r  s  o 


a.  =  a  (0.25  +  5.5  H  /L  )  (i=l,2) 

D  .  O  .  O  O 


* 


where : 


r  \  cos  a. 


“•s  (v^)' 


o  1.418 


Note  also  that  in  the  above  equations,  the  subscripts  1  and  2  refer  to 


the  Xj  and  contour  lines,  respectively. 


Depth  of  Closure 


Following  the  approach  of  Kraus  and  Harikai  (1983),  the  depth  of 
closure  model  developed  by  Hallermeier  (1978)  was  used  to  describe  the 
seaward  boundary  of  the  zone  where  significant  longshore  and  intense 
on-offshore  transport  takes  place.  Hallermeier ' s  model  expresses  the 


critical  depth,  as 


D  =  2.28  H  -  68.5  (H  ^/gT^) 

c  s  s  ^  s 


where  and  are  expressed  in  meters,  g  in  m/sec  ,  and  in  seconds. 


Boundary  and  Initial  Conditions 


The  initial  conditions  for  the  model  are  completely  arbitrary, 
however  for  simplicity,  the  shoreline  is  assumed  to  be  initially  parallel 


to  the  y  axis,  with 
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and 


^2  ”  ^  Dl/tan  (46) 

where:  tan  =  the  initial  beach  slope 

The  boundary  conditions  for  the  model  are  specified  at  both  the 
groin  and  the  far  upcoast  and  downcoast  positions.  At  the  latter 
positions,  it  is  assumed  that 


9  X  . 
_ 1 
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0  i=l,2 


(47) 


At  the  groin,  the  boundary  condition  is  specified  in  terms  of  the  rate 
of  longshore  transport.  It  is  assumed  that  sand  is  allowed  to  bypass 
the  groin  in  proportion  to  r|  ,  which  is  defined  as  follows.  Letting  the 
dimensionless  groin  length,  x*  =  h^/h^,  where  h^  is  the  depth  of  water 
at  the  end  of  the  groin,  then  the  coefficient  r)g  can  be  evaluated  using 
Equation  42.  Having  specified  the  boundary  condition  at  the  groin 
is  specified  as 


If  n  <  = 

g 


Hj  then: 


n,  =  n,  -  Hg  .nd 


If  n,  <  Og  < 


then; 


(48) 
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If  n. 


then : 


I  t 

=  0  and  =  0 

where:  the  superscript  '  =  value  of  ^2  groin 


Finite  Difference  Solution 


Following  the  approach  of  Perlin  and  Dean  (1978),  an  explicit 
solution  utilizing  a  time-marching  space  and  time-staggered  step 
procedure  was  used.  This  procedure  involves  holding  the  shoreline 
orientation  fixed  for  one-time  step  (from  nt  to  (n+l)t)  and  calculating 
the  longshore  and  on-offshore  sand  transports.  The  sand  transport  is 
then  held  fixed  over  a  time  step  (from  (n+l/2)t  to  (n+3/2)t)  and  the 
shoreline  changes  are  determined.  The  finite  difference  equations  are 


>r‘ 

1 


n+1  n+1 

Hi  K. 


n+1/2  n+1/2 

sin  a.  cos  a. 


(49) 


1 . 
1 


1 . 
1 


n+1 

^2 

i 


1  .,n+l  '  'i  i 


n+.  „ 

^12  •'i 


n+1/2  n+1/2 

8  /S  -1\  N  “b„  "b„ 

V  s  /  o  2^  2^ 


(50) 


<i  (k  h",  tan 

X  \  R  s  / 1 


X  . 
1 


p  g 


(51) 


n+1  n+1 

-Q  +Ayq 


n+3/2 

x 


n+1/2 

X, 


At _ 

/'D_+  D.)  Av  '^l. 


(52) 


n+3/2 

- 

At 

n+1  n+1  n+1 

(53) 

2  . 

"  2. 

Do 

1 

1 

2 

1+1  1  1 

n+1/2 
b  . 

Ji 

n+1/2 

=  a 

o  . 

Ji 

^0.25  +  : 

i.5  j=l,2 

(54) 

,  a.l/2  ^ 

n+1/2  -1  V  ^i  ^i-l  / 

j  =  l,2 

(55) 

=  or  -  tan  - r - 
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where  the  superscripts  denote  the  time  level  at  which  the  variable  is 
evaluated  and  the  subscripts  denote  the  space  step  at  which  the  variable 
is  evaluated. 


Example  Solutions 

Two  examples  are  presented  which  illustrate  the  shoreline  evolution 
around  a  newly  constructed  groin  superimposed  on  an  overall  changing 
beach  profile.  In  the  first  example,  the  initial  beach  slope  is  steeper 
than  the  equilibrium  beach  slope.  In  the  second  example,  the  initial 
beach  slope  is  flatter  than  the  equilibrium  beach  slope.  Input  conditions 
for  the  two  examples  are  as  follows: 

Example  1: 

=  2.12  meters;  T  =  8  seconds;  or  =  -25  degrees 

tan  p.  =  0.02;  D,  =  1.0  meters;  D,  =  1.5  meters 

1  b  1 

W  =  0.04  m/sec;  Ay  =  100  meters;  At  =  3  hours 

X  =  100  meters 

groin 

Example  2: 

H 

s 


0.71  meter;  T 


8  seconds;  a 


-25  degrees 


0.75  meters 


tan  =  0.15;  =  1.0  meters; 


W  =  0.04  m/sec;  Ax  =  100  meters;  At  =  12  hours 


X  .  -  SO  meters 

groin 


Figures  16  and  17  show  the  time  evolution  of  the  nearshore  contours  for 
Examples  1  and  2  respectively.  Contour  positions  are  plotted  for  times 
of  0,  1,  4,  16,  and  64  days.  Although  the  solution  is  shown  for  both 
sides  of  the  groin,  the  contour  line  positions  on  the  immediate  down-draft 
side  of  the  groin  are  not  realistic  because  wave  diffraction  has  been 
neglected . 

In  the  immediate  vicinity  of  the  groin.  Figures  16  and  17  show 
similar  patterns  of  shoreline  behavior,  with  a  sand  fillet  forming  on 
the  upcoast  side  of  the  groin,  and  erosion  occurring  on  the  downcoast 
side.  In  addition,  the  offshore  contours  show  a  steepening  of  the  beach 
on  the  upcoast  side  of  the  groin,  and  a  flattening  of  the  beach  on  the 
downcoast  side.  This  behavior  is  in  accordance  with  field  observations. 

Superimposed  on  the  above  general  trends  is  an  overall  flattening 
of  the  entire  beach  in  Figure  16  caused  by  a  general  offshore  sand 
transport.  In  Figure  17,  the  beach  is  slowly  steepening  with  a  general 
onshore  sand  transport.  The  effect  of  these  global  changes  is  to 
exaggerate  the  effects  of  the  local  accretion  and  erosion  around  the 
groin  in  Example  1  and  diminish  the  effects  in  Example  2. 


DISCUSSION 


The  primary  advantages  of  the  longshore  and  on-offshore  transport 
models  developed  in  this  study  are  that  they  share  a  common  energetics 
development  and  they  are  simple  to  apply.  The  equation  for  the  wave 
power  coefficient,  K,  involves  parameters  which  are  readily  available 
from  the  usual  longshore  transport  calculations.  Separate  accounting  of 
the  bedload  and  suspended  load  transports  provides  for  a  greater  degree 
of  flexibility  in  the  wave  power  equation,  allowing  the  equation  to  be 
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used  over  a  wider  range  of  wave  heights  and  sediment  sizes.  The 
on-offshore  transport  model  is  extremely  easy  to  use,  requiring  only  the 
incident  wave  height,  the  sediment  fall  velocity,  and  an  average  drag 
coefficient  for  the  bed.  When  combined  with  a  simple  two-line  shoreline 
evolution  model,  the  above  models  can  be  used  to  predict  beach  profile 
and/or  shoreline  evolution  around  a  groin. 

The  above  models  have  several  general  and  specific  limitations. 

The  general  limitations  are  primarily  related  to  gross  simplifications 
made  in  deriving  the  energetics  transport  equations.  First,  the 
longshore  and  on-offshore  sediment  transport  models  are  vertically 
integrated,  thus  ignoring  the  vertical  structure  of  the  sediment 
concentration  and  velocity  fields.  Recent  field  measurements  (Jaffe,  et 
al.,  1984)  suggest  that  the  vertical  distributions  of  these  parameters 
are  important. 

Another  general  limitation  of  the  energetics  models  is  that  the 
rate  of  energy  dissipation  (to  which  the  transport  rate  is  proportional) 
is  explicitly  related  to  the  bottom  boundary  layer  (i.e.,  the  drag  on 
the  bottom).  Contributions  to  the  local  rate  of  energy  dissipation  from 
the  breaking  wave  at  the  surface  are  neglected,  in  spite  of  the  fact 
that  the  latter  has  been  shown  to  be  dominant  (Thornton  and  Guza ,  1983). 
An  implicit  assumption  used  in  the  present  model  development  was  that 
the  bottom  boundary  layer  controls  the  upward  flux  (entrainment)  of 
sediment.  This  assumption,  however,  is  only  conjecture.  Possible 
improvements  in  the  model  might  be  to  incorporate  a  variable  drag 
coefficient  to  account  for  some  of  the  effects  of  breaking  waves; 
alternatively,  the  rate  of  energy  dissipation  could  be  directly 
associated  with  the  local  gradient  in  the  energy  flux  (e.g.,  Thornton, 
1973) . 

Ignoring  sediment  threshold  of  motion  effects  is  another  general 
limitation  of  the  models.  In  the  case  of  longshore  transport  of  sand¬ 
sized  sediments,  under  most  field  conditions  threshold  effects  are 
negligible.  However,  threshold  effects  are  probably  significant  in 
on-offshore  transport  under  most  field  conditions.  Similarly,  they  are 
probably  important  for  longshore  transport  under  laboratory  conditions. 
The  statistical  method  developed  by  Seymour  (unpublished  manuscript, 
1984)  may  be  one  method  of  dealing  with  these  effects. 


Specific  model  shortcomings  can  be  summarized  as  follows.  The 

longshore  transport  distribution  was  derived  using  common  spilling  wave 
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assumptions  (e.g.,  the  orbital  velocity  magnitude  decreases  as  h  ). 

Field  measurements  on  dissipative  beaches,  however,  have  shown  that  the 
total  velocity  variance,  and  thus  the  orbital  velocity  magnitude,  is 
nearly  constant  across  the  surfzone.  The  longshore  current  model  used 
in  the  present  study  is  also  largely  unverified.  Beach  profile  shapes 
have  been  shown  to  have  a  signficant  effect  on  the  distribution  of  the 
longshore  current  (Hudspeth  and  McDougal,  1984).  Moreover,  the  proper 
formulation  for  the  lateral  mixing  of  momentum  is  still  an  unanswered 
question.  As  a  result,  the  distribution  of  the  longshore  transport 
rate,  as  specified  in  the  present  model,  should  be  viewed  as  approximate. 
The  spatially  integrated  equation  for  the  wave  power  coefficient,  K,  is 
probably  more  accurate,  because  uncertainties  in  the  transport  distribution 
have  been  hidden  in  the  spatial  integration. 

The  present  on-offshore  transport  model  is  only  a  first  step  in 
describing  the  on-offshore  transport.  Although  it  is  dynamically  based, 
the  empirically-derived  regression  equations  for  the  wave  velocity 
moments  are  restricted  to  a  narrow  range  of  conditions.  In  fact,  the 
present  model  is  unable  to  predict  onshore  sand  motion  for  beach  slopes 
in  excess  of  approximately  0.035.  Eventually,  when  theoretical  or 
numerical  models  are  developed  to  predict  the  wave  velocity  moments  as  a 
function  of  the  incident  wave  field  and  the  local  water  depth,  it  will 
be  possible  to  predict  the  time-evolution  of  the  entire  beach  profile. 

This  would  allow  a  direct  comparison  of  the  present  dynamically  based 
on-offshore  transport  model  with  Swart's  (1974)  empirically  based  beach 
profile  model. 

The  line  model  methodology  used  in  the  present  study  is  relatively 
well  established.  The  contributions  of  this  study  relate  to  the 
incorporation  of  improved  longshore  and  on-offshore  sediment  transport 
models.  A  further  improvement  of  the  present  line  model  would  be  to 
include  wave  diffraction  effects,  using  the  techniques  developed  by 
LeMehaute  and  Soldate  (1980)  or  Kraus  and  Harikai  (1983).  Extension  of 
the  model  to  three  or  more  lines  is  not  feasible  at  this  time  due  to  the 
limitations  in  the  on-offshore  transport  model. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  conclusions  were  reached  as  a  result  of  this  study. 

1.  An  energetics-based  longshore  sediment  transport  model  was  developed 
which  resulted  in  a  modification  to  the  wave  power  coefficient.  Instead 
of  being  const. int,  K  was  found  to  be  a  function  of  the  breaker  angle, 
the  beach  slope,  and  the  ratio  of  the  orbital  velocity  magnitude  at  the 
breakpoint  divided  by  the  fall  velocity  of  the  sediment.  These  modifi¬ 
cations  were  found  to  extend  the  range  of  applicability  of  the  wave 
power  equation  to  both  laboratory  and  field  conditions.  Although  the 
model  describes  the  distribution  of  the  longshore  transport  rate  across 
the  surfzone,  this  distribution  should  be  applied  with  caution.  Recent 
measurements  of  surfzone  wave  heights  and  currents  suggest  that  a  number 
of  common  assumptions  used  in  the  development  of  the  present  model,  and 
most  longshore  current  models,  are  not  accurate. 

2.  An  energetics-based  on-offshore  sediment  transport  model  was  developed 
which  predicts  the  spatially  averaged  surfzone  on-offshore  transport 

rate  as  a  function  of  the  incident  wave  height,  the  average  beach  slope, 
and  the  sediment  fall  velocity.  While  the  model  remains  unverified,  it 
exhibits  a  qualitative  behavior  which  corresponds  to  field  observations 
of  beach  profiles.  Additional  work  is  needed  to  verify  the  existing 
model  and  to  extend  the  model  to  predict  the  cross-shore  distribution  of 
the  on-offshore  sediment  transport  rate. 

3.  The  above  two  sediment  transport  models  were  used  to  develop  a 
two-line  model  which  describes  the  evolution  of  the  shoreline  and  beach 
slope  as  a  function  of  the  incident  wave  conditions.  The  model  describes 
both  the  average  beach  profile  changes  (i.e.,  the  shoreline  position  and 
the  average  beach  slope)  along  an  open  beach  as  well  as  the  shoreline 
changes  associated  with  a  single  groin.  Additional  work  is  needed  to 
incorporate  wave  diffraction  effects. 


4.  There  is  a  limit  to  the  extent  that  vertically  integrated 
energetics-based  sediment  transport  models  can  be  used  to  describe 
surfzone  sediment  transport  processes.  The  present  models  can  be 
modified  by  including  sediment  threshold  of  motion  effects,  a  variable 
drag-coefficient,  and  a  different  formulation  for  the  rate  of  energy 
dissipation  in  the  surfzone.  These  modifications  will  probably  result 
in  an  improvement  to  the  models,  especially  in  the  case  of  the  on-offshore 
sediment  transport  model.  Most  likely,  however,  significant  progress  in 
developing  improved  sediment  transport  models  will  only  come  through  an 
increase  in  our  understanding  of  the  vertical  structure  of  the  sediment 
concentration  and  velocity  fields. 

Based  on  these  findings,  it  is  recommended  that  future  work  focus 
on  measuring  the  vertical  structure  of  the  sediment  concentration  and 
fluid  velocity  fields.  A  parallel  effort  should  focus  on  developing 
improved  models  for  wave  bore  dynamics,  surfzone  boundary  layer  development, 
bedload,  and  suspended  load  mechanics.  Improved  models  for  the  longshore 
and  on-offshore  currents  in  the  surfzone  are  also  needed. 
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Table  1 .  Longshore  Transport  Data 
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Silver  Strand 

3 

2.0 

72 

4 

0.66 

Kraus  et  al.  (1983) 
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0.31 

Hirono 
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0.11 

Shimokita 
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6.0 

0.61 

Moore  &  Cole  (I960) 
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^  *  (ECn)j,  sinoj,  cosa^  .  erg/sec  cm 

Fiii’irt'  I.  The  wave  pov^er  eauation  describes  the  field 
measurements  of  lonpshore  transport  re]ativelv  well, 
however,  it  over-predicts  the  laboratory  measurements 
(from  P.D.  Komar  and  D.L.  Inman  (1970),  "Lonpshore 
sand  transport  on  beaches,"  Journal  of  Geophysical 
Research,  vol  75,  no.  30,  1970,  Fig.  5,  pg  5922). 
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Figure  2.  Progressive  accretion  and  erosion  near  a 
groin,  based  on  a  one-line  numerical  model.  The 
solid  lines  represent  the  solution  with  diffrac¬ 
tion,  the  dashed  lines  without  (from  W.T.  Bakker 
K.  Breteler,  A,  Roos,  1970,  "The  dynamics  of  a 
coast  with  a  groyne  svstem,"  ASCE  Proc.  12th 
Coastal  Eng,  Conf.,  Fig.  5,  pg  1004,  vol  2). 


Figure  3.  Definition  sketch  for  two-line 
shoreline  evolution  model. 


Lon^hore  Transport 
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Figure  4.  Evolution  of  the  shoreline  and  offshore 
beach  limit  near  a  groin  (from  W.T. 

Bakker,  1968,  "The  dynamics  of  a  coast 
with  a  groin  system,"  ASCE  Proc.  11th 
Coastal  Eng.  Conf,,  Fig.  8,  pg  497,  vol  1) 


Daan  et.al 


Figure  6.  Observed  values  of  the  wave  power  coefficient  versu 
estimated  values  using  Equation  33.  The  solid  line 
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Figure  7.  Significant  wave  height  and  wave  period  characteri 
for  the  Torrey  Pines  Beach  and  Leadbetter  Beach  da 
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Appendix  A 


A  DERIVATION  OF  AN  EXPRESSION  FOR  THE  BEACH  SLOPE 
EXPRESSED  AS  A  FUNCTION  OF  TIME,  BASED  ON  THE 
BARKER/ SWART  ON-OFFSHORE  TRANSPORT  MODEL 


Assuming  zero  divergence  of  the  longshore  sediment  transport  rate,  the 
two-line  sand  conservation  equations  (Equations  6  and  7)  become 


(A-1) 


a  X, 

)  - - 

^2  at 


=  <1. 


(A-2) 


where  q  ,  the  on-offshore  sediment  transport  rate  is  defined  as 

q^  =  Coff  ”  *2  (A“3) 

and  w  =  the  equilibrium  value  of  (x^  -  x^).  Defining 

z  =  Xj  -  x^  +  w  (A-4) 

and  taking  the  derivative  of  z  with  respect  to  t,  the  following  result 
is  obtained 


_1 _ 

Coff 


_ 


(A-5) 


A-1 


-i'.- 


where  use  has  been  made  of  Equations  A-1  and  A-2.  Rearranging 
Equation  A-4,  and  integrating,  the  following  equation  for  the  on-offshore 
transport  rate  is  obtained 


(-A6) 


where  C  is  an  integration  constant,  and  B 


Coff  Coff 
°2  ■ 


Applying  the  following  initial  condition  (i.e.,  at  time  t  =  0,) 


q  =  Coff 

X 


tan  6 . 
1 


2  D 


1 


tan  p 


(A-7) 


where  tan  p  is  the  initial  beach  slope,  and  tan  p  is  the  equilibrium 
beach  slope,  the  constant  C  becomes 


C  =  Coff 


-2  D 


1 


2  D 


1 


tan  p.  tan  p 


(A-8) 


Substituting  Equation  A-7  into  Equation  A-5,  and  rearranging,  the  desired 
equation  for  the  beach  slope,  tan  p,  is  obtained 


Appendix  B 


A  DERIVATION  OF  AN  EQUATION  FOR  Hj 

In  a  two-line  model  is  it  necessary  to  split  the  longshore  transport 
rate  into  two  parts,  that  which  occurs  shoreward  of  the  Dj  depth  contour, 
and  that  which  occurs  seaward  of  the  depth  contour.  The  parameter, 

Hj,  is  defined  as  the  fraction  of  the  total  longshore  transport  rate 
occurring  between  the  shoreline  and  the  D-  contour  line.  The  normalized 
X  coordinate  position  of  this  contour  line  x^  is  defined  as 


* 


where : 


In  this 


h^  =  water  depth  at  the  break  point 
way,  r).j^  can  be  expressed  as 


?1  " 
f  <i  >dx 
•'o  y 


/  <i  >dx 

y 

o 


(B-1) 


(B-2) 


Recalling  that  Equation  22,  which  defines  the  cross-shore  distribution 
of  the  longshore  sediment  transport  rate,  assuming  that  the 

transport  is  primarily  by  suspension,  then  following  a  lengthy 
integration,  the  following  expression  is  obtained 

* 

If  X  <  1  then 


[f  Xj  >  1  then 


"i  =  ’  *57^-57^  I, 


'  -  s  \  s  - 


Finally,  if  P  =  0.2  then 

.4.11 


,3.5 


fl  j  =  -1.69  Xj  +  2.34  Xj 


.-3.61 


1  -  0.35  X, 
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ADMINSUPLi  PWO,  Bahrain 

BUREAU  OF  RECLAMATION  Code  1.512  (C.  Sclander)  Denver  CO 

CBC  Code  10.  Davisville.  RI;  Code  L55,  Port  Hueneme,  CA;  Code  156.  Port  Hueneme.  CA;  Code  156F,  Port 
Huenemc.  Ca;  Code  1571,  Port  Hueneme.  CA;  Dir,  CESO,  Port  Hueneme.  CA;  Library.  Davisville,  RI; 
PWO.  Gulfport,  MS;  Tech  Library.  Gulfport,  MS 
CBU  401.  OICC.  Great  Lakes,  IL;  411,  OIC,  Norfolk,  VA 
CINCLANTFLT  Civil  Ei.j,  Supp  Plans  Offr,  Norfolk,  VA 

CNO  Code  NOP-964,  Washington  DC;  Code  OP  323.  Washington  DC;  Code  OP  414.  Washington  DC;  Code  OP 
424,  Washington  DC;  Code  OP  97  Washington  DC;  Code  OP  97  Washington.  DC;  Code  OP  987  Washington 
DC,  Code  OP-413  Wash.  DC;  Code  OP-987J.  Washington.  DC;  Code  OP323  Washington  DC;  Code 
OPNAV  09B24  (H);  Code  OPNAV  22,  Wash  DC;  Code  OPNAV  23.  Wash  DC;  OP-41  IF.  Wash  DC 
COMCBLANT  Code  S.3T,  Norfolk,  VA 
COMCBPAC  Diego  Garcia  Pro)  Offr.  Pearl  Harbor.  HI 
COMFAIRMED  SCE  (Code  N55),  Naples,  Italy 
COMFEWSCi  DET  Security  Officer  (R.  Seidman).  Washington,  DC 
C  OMFLEACT,  OKINAWA  PWO.  Kadena,  Japan 
COMNAVACT  PWO.  London.  England 
COMNAVAIRSYSCOM  Code  41712.  Washington,  DC 

COMNAVBEACHGRU  ONE,  CO.  San  Diego,  CA;  TWO,  CO.  Norfolk.  VA 
CfJMNAVFORCARIB  Code  N42.  Roosevelt  Rds.  PR 
COMNAVFORKOREA  ENJ-P&O.  Yongsan 
COMNAVLOGPAC  Code  4318,  Pearl  Harbor.  HI 
COMNAVMARIANAS  CO.  Guam 
COMNAVMEDCOM  Security  Officer.  Washington.  DC 
COMNAVRESFOR  Code  08.  New  Orleans.  LA 
COMNAVSUPPFORANTARCTICA  DET.  PWO,  Christchurch,  NZ 
COMNAVSURFLANT  CO.  Norfolk,  VA;  Code  N42A  Norfolk.  VA 
COMNAVSURFPAC  Code  N-4,  San  Diego.  CA 
COMOCEANSYSPAC  SCE,  Pearl  Harbor,  HI 
COMSC  Washington  DC 

COMSPAWARSYSCOM  Code  PME  124-61,  Washington,  DC;  PME  124-612,  Washington,  DC 
COMSUBDEVGRUONE  Operations  Offr,  San  Diego,  CA 
COMSURFWARDEVGRU  CO,  Norfolk,  VA 
COMUSNAVCENT  Code  N42,  Pearl  Harbor,  HI 


NAVOCEANC'OMC'EN  CO,  Guam,  Mariana  Islands;  Code  EES,  Guam.  Mariana  Islands 
NAVRESCEN  PE-PLS,  Tampa,  FL 

COMOPTEVFOR  CMDR,  Norfolk.  VA,  Code  H)5.  San  Diego.  CA 

DEFFUELSUPPCEN  DFSC-OWE  (Term  Engrng)  Alexandria,  VA;  DFSC-OWE.  Alexandria  VA 

DI.A  DB-6E1.  Washington.  DC;  DB-6E2,  Washington,  DC 

DLSIE  Army  Logistics  Mgt  Center.  Fort  Lee.  VA 

D.NA  SITI'TL.  Washington,  DC 

DOE  Wind'Oeean  Tech  Div.  Washington.  DC 

DTIC  Alexandria,  VA 

DTNSRDC  Code  IP.S,  Annapolis.  MD;  Code  119.  Annapolis,  MD;  Code  1250.  Annapolis,  MD;  Code  L561. 
Bethesda.  MD;  Code  1.568,  Annapolis,  MD;  Code  1706,  Bethesda  MD;  Code  172,  Bethesda.  MD;  Code  284, 
Annapolis.  MD;  Code  2842.  Annapolis,  MD;  Code  4111  (R.  Gierich).  Bethesda  MD;  Code  4120,  Annapolis, 
MD 

EODGRU  ONE  DET,  CO.  Point  Mugu,  CA 
FAA  (Fowler)  Code  APM-740,  Wash,  DC 
FMFLANT  CEC  Offr.  Norfolk  VA 

FMFPAC  FEO,  Camp  HM  Smith,  HI;  G5  (SCIAD),  Camp  HM  Smith,  HI 
FOREST  SERVICE  Engrg  Staff,  Washington.  DC 
GIDEP  OIC,  Corona,  CA 

GSA  Assist  Comm  Des  &  Cnst  (FAIA)  D  R  Dibner  Washington,  DC 
IRE-FITD  Input  Proc  Dir  (R.  Danford).  Eagan,  MN 
KWAJALEIN  MISRAN  BMDSC-RKL-C 

LIBRARY  OF  CONGRESS  Washington,  DC  (Sciences  &  Tech  Div) 

MARCORDIST  12.  Code  4,  San  Francisco,  CA 

MARCORPS  FIRST  FSSG,  Engr  Supp  Offr,  Camp  Pendleton.  CA 

MARCORPS  AIR/GND  COMBAT  CTR  ACOS  Fac  Engr,  Okinawa,  Japan 

MARINE  CORPS  BASE  ACOS  Fac  engr,  Okinawa;  Dir,  Maim  Control.  PWD.  Okinawa.  Japan;  PWO,  Camp 
Lejeune,  NC;  PWO,  Camp  Pendleton  CA 

MARINE  CORPS  HOTRS  Code  LFF-2,  Washington  DC;  Code  LM-2,  Washington,  DC 
MCAS  Facil.  Engr,  Div.  Cherry  Point  NC 
MCAF  Code  C144.  Ouantico.  VA 

MCAS  Facs  Maint  Dept  -  Operations  Div,  Cherry  Point;  PWO,  Kaneohe  Bay,  HI;  PWO,  Yuma  AZ 
MCDEC  M  <&  L  Div  Ouantico,  VA;  NSAP  REP,  Ouantico  VA 
MCRD  SCE  San  Diego  CA 
NAF  PWO.  Atsugi  Japan 

NAS  Code  16.1.  Kcflavik,  Iceland;  Code  182,  Bermuda;  Code  187,  Jacksonville  FL;  Code  22.  Patuxent  River. 
MD;  Code  81.  Patuxent  River.  MD;  Code  SEN,  Patuxent  River,  MD;  Dir,  Engrg  Div,  Millington,  TN; 
Director,  Engrg,  Div;  Engrg  Dir,  PWD,  Corpus  Christi,  TX;  P&E  (Code  1821H).  Miramar,  San  Diego,  CA; 
Code  72E,  Willow  Grove,  PA;  Lead  CPO,  PWD,  Self  Help  Div,  Beeville,  TX;  PWD  Maint  Div,  New 
Orleans,  LA;  Code  72E.  Willow  Grove,  PA;  PWO  New  Orleans.  LA;  PWO,  Glenview  IL;  PWO,  Keflavik, 
Iceland;  PWO.  Key  West,  FL;  PWO,  Milton,  FL;  PWO,  Sigonella.  Sicily;  PWO,  South  Weymouth,  MA; 
PWO,  Willow  Grove,  PA;  SCE  Norfolk,  VA;  SCE,  Barbers  Point,  HI;  SCE,  Cubi  Point,  RP 
NATL  BUREAU  OF  STANDARDS  R  Chung,  Gaithersburg,  MD 
NATL  RESEARCH  COUNCIL  Naval  Studies  Board,  Washington  DC 
NAVAIRDEVCEN  Code  81.1.  Warminster  PA;  Code  8323,  Warminster,  PA 
NAVAIRENGCEN  PWO.  Lakehurst.  NJ 

NAVAIREWORKFAC  Code  100,  Cherry  Point,  NC;  Code  640,  Pensacola  FL;  Equip  Engr  Div  (Code  61000), 
Pen,sacola,  FA 

NAVAIRTESTCEN  PWO.  Patuxent  River,  MD 

NAVCAMS  PWO.  Norfolk  VA;  SCE  (Code  N-7),  Naples,  Italy;  SCE,  Wahiawa  HI;  SCE.  Wahiawa,  HI; 
Security  Offr,  Wahiawa,  HI 

NAVCHAPGRU  CO  Williamsburg  VA;  Operations  Officer.  Code  30  Williamsburg,  VA 

NAVCOASTSYSCEN  CO,  Panama  City.  FL;  Code  2230  (J.  Quirk)  Panama  City,  FL;  Code  715  (J.  Mittleman) 
Panama  City,  FL;  Code  772  (C.B.  Koesy),  Panama  City,  FL;  PWO,  Panama  City,  FL;  Tech  Library, 

Panama  City,  FL 

NAVCOMMSTA  Code  401,  Nea  Makri.  Greece;  Dir,  Maint  Control.  PWD,  Diego  Garcia;  PWO,  Exmouth, 
Australia 

NAVCONSTRACEN  Code  0OU15,  Port  Hueneme  CA;  Code  B-l,  Port  Hueneme,  CA;  Curriculum/Instr.  Stds 
Offr,  Gulfport  MS 

NAVEDTRAPRODEVCEN  Tech  Library.  Pensacola.  FL 
NAVEODTECHCEN  Tech  Library.  Indian  Head,  MD 

NAVFAC  Maint  &  Stores  Offr.  Bermuda;  PWO.  Centerville  Bch,  Ferndale  CA 

NAVFACENGCOM  CO,  Alexandria,  VA;  Code  03,  Alexandria,  VA;  Code  03T  (Essoglou),  Alexandria,  VA; 
Code  04A1.  Alexandria,  VA;  Code  04B3,  Alexandria,  VA;  Code  04M,  Alexandria,  VA;  Code  04T1B 
(Bloom).  Alexandria,  VA;  Code  04T4,  Alexandria.  VA;  Code  04X5,  Alexandria.  VA;  Code  06,  Alexandria 
VA;  Code  07A  (Herrmann),  Alexandria,  VA;  Code  09M124  (Tech  Lib).  Alexandria,  VA;  Code  100, 
Alexandria,  VA;  Code  1113,  Alexandria,  VA 


NAVFAC'FNGC'OM  ■  CUES  DIV  Ciidc  4115,  Washington,  DC;  Code  407  (D  Schccscic)  Washington,  DC;  Code 
FPO-IC  Washington  DC;  FPO-1  Washington,  DC;  Library,  Washington,  D  C 
NAVFACENCiCOM  -  LANE  DIV  Br  Otc,  Dir,  N;iples,  Italy;  Code  1117,  Norfolk,  VA;  Code  40.4,  Norfolk, 
VA;  Code  405.  Norfolk.  VA;  Library.  Norfolk.  VA 
NAVFACENCiCOM  -  NORTH  DIV.  C(7.  Philadelphia.  PA;  Code  04AL.  Philadelphia  PA;  Code  (WP. 

Philadelphia.  PA;  Code  11,  Philadelphia,  PA;  ROICC.  Contraets.  C  rane  IN 
NAVFACENCiCOM  -  PAC  DIV.  Code  (NP.  Pearl  Harbor.  HI;  Code  7011  Pearl  Harbor.  HI;  Code  407, 
RDT&E.  Pearl  Harbor.  HI,  Library.  Pearl  Harbor,  HI 
NAVFACENCiCOM  -  SOUTH  DIV.  Code  1112,  Charleston.  SC;  Code  406  Charleston.  SC;  Library, 

Charleston.  SC 

NAVF’AC  EN'CiC OM  -  WEST  DIV.  (WP  20,  ,San  Bruno,  CA;  Code  (I4B.  San  Bruno.  CA:  Library.  San  Bruno. 
C.A;  RDl&E  LnO.  San  Bruno,  CA 

NA\  FAC  ENCiCOM  CONTRACTS  DOICC.  Diego  Garcia;  OICC  .  Guam;  OICC  .  Rota  Spain;  OICC  RCdICC  . 
Norfolk.  VA;  ROICC  (Code  405).  Portsmouth,  VA;  ROICC.  Corpus  Chrisfi.  IX;  ROICC.  KeBavik. 

Iceland;  ROICC.  Key  West,  FL.  ROICC.  Point  Mugu.  CA;  ROICC,  Iwcntynine  Plams.  CA;  ROICC  OIC  C  . 
SPA,  Norfolk.  VA;  SW  Pac,  Dir,  Engr  Div.  Mama.  RP;  SW  Pac.  OICC,  Manila,  RP;  1,'ident,  OICC.  St 
Marvs.  GA 

NAS'HOSP  CO.  Long  Beach.  CA.  PWO.  Guam.  Mariana  Islands;  SCE.  Pensacola  FL,  SCE.  Yokosuka.  Japan 
NAVMAG  SCE.  Subic  Bay  RP 

NAVMEDCCJM  SEREG.  Head.  Fac  Mgmt  Dept,  Jacksonville.  FL 

N.AVOCEANO  Code  .44,47  (J  DePalma).  Bay  St.  Louis  MS;  Code  627(1  (M.  Paige),  Bay  St.  Louis.  MS;  Library 
Bay  St.  Louis.  MS 

NAVOCEA.NSYSCEN  Code  52(14  (J.  Stachiw),  San  Diego.  CA;  Code  541  (Bachman)  San  Diego.  CA;  Code  90 
(Talkington),  San  Diego,  CA;  Code  944  (H.C  Wheeler).  San  Diego.  CA;  C'ode  964  (Tech  Library),  San 
Diego,  CA;  Code  9642B  (Bay.side  Library).  San  Diego.  CA;  DET.  R  Yumori.  Kailua.  HE.  DET.  Tech  Lib. 
Kailua.  HI 

NAVPETOFF  Code  8D1(I7.  Alexandria.  VA 
NAVPETRES  Director,  Washington  DC 

NAVPGSC'OL  C.  Morers,  Monterey,  CA;  Code  1424.  Library.  Monterey.  CA;  Code  61WL  (O.  Wilson), 
Monterey,  CA;  C'ode  6K  (C.S.  Wu).  Monterey.  CA;  E.  Thornton.  Monterey.  C'A 
NAVPHIBASE  Harbor  Clearance  Unit  Two,  Norfolk,  VA;  PWO  Norfolk,  VA;  SCE.  San  Diego.  CA 
NAVRESREDCOM  Commander  (Code  072),  San  Francisco.  CA 
NAVSCOLCECOFF  C.45  Port  Hueneme,  CA 
NAVSCSCOL  PWO.  Athens  GA 

NAVSEASYSCOM  Code  (143,  Washington  DC;  Code  0,5G14.  Washington,  DC;  Code  05R12,  Prog  Mgr 
Washington.  DC;  Code  ()6H4,  Washington,  DC;  Code  CL42.  Washington.  DC;  Code  PMS-.496.421 1  (J. 

Rekas)  Washington.  DC 

NAVSECGRUACT  PWO  (Code  .4(15).  Winter  Harbor,  ME;  PWO  (Code  40).  Edzell.  Scotland 
NAVSECGRUCOM  Code  G44.  Washington,  DC 

NAVSHIPYD  Carr  Inlet  Acou.slic  Range,  Bremerton.  WA;  Code  202.5  (Library).  Bremerton.  WA;  Code  280. 
Mare  Is..  Vallejo.  CA:  Code  280.28  (Goodwin),  Vallejo.  CA;  Code  480.  Portsmouth.  VA;  Code  410.  Mare 
Is..  Vallejo  CA;  Code  440,  Bremerton,  WA;  Code  440.  Bremerton,  WA;  Code  440,  Norfolk.  VA;  Code  440. 
Portsmouth.  NH;  Code  440.4,  Bremerton,  WA;  Dir,  Maint  Control.  PWD,  Long  Beach.  CA;  PWO, 
Charleston.  SC;  PWO,  Mare  Island,  Vallejo,  CA 

NAVSTA  A.  Sugihara,  Pearl  Harbor.  HI;  CO.  Long  Beach,  CA;  CO.  Roosevelt  Roads.  PR;  Dir  Mech  Engr 
47WC94  Norfolk,  VA;  Dir,  Engr  Div,  PWD  (Code  I82(X)).  Maypori,  FL;  Engrg  Dir.  Rota.  Spain;  PWO. 
Mayport.  FL;  SCE,  Subic  Bay.  RP;  Util  Engrg  Offr.  Rota,  Spain 
NAVSUPPFAC  Dir,  Maint  Control  Div,  PWD,  Thurmont.  MD 
NAVSUPPO  Security  Offr.  La  Maddalena,  Italy 

NAVSURFWPNCEN  Code  E211  (C.  Rouse),  Dahlgren.  VA;  DET.  Security  Offr,  Silver  Springs,  MD;  G-.S2 
(Duncan)  Dahlgren.  VA;  PWO,  Dahlgren.  VA 
NAVTECHTRACEN  SCE,  Pensacola  FL 
NAVWARCOL  Fac  Coord  (Code  24).  Newport,  Rl 

NAVWPNCEN  Code  2646,  China  Lake.  CA;  DROICC  (Code  702).  China  Lake.  CA 

NAVWPNSTA  Dir.  Maint  Control.  PWD,  Concord,  CA;  Engrg  Div,  PWD.  Yorktown.  VA;  PWO.  Charleston, 
SC;  PWO,  Code  096.  Colts  Neck.  NJ;  PWO,  Seal  Beach,  CA 
NAVWPNSTA  PWO,  Yorktown.  VA 
NAVWPNSUPPCEN  Code  09  Crane  IN 
NETC  PWO.  Newport.  Rl 
COMEODGRU  OIC,  Norfolk  VA 
NCR  20.  CO.  Gulfport,  MS;  20.  Code  R70 

NMCB  74.  CO;  FIVE.  Operations  Dept;  Forty,  CO;  THREE,  Operations  Off, 

NOAA  (Mr.  Joseph  Vadus)  Rockville,  MD;  Library,  Rockville.  MD 

NORDA  Code  410.  Bay  St.  Louis,  MS:  Head,  Geotech  Br  (Code  .463),  Bay  St.  Louis,  MS;  Ocean  Rsch  Off 
(Code  440).  Bay  St.  Louis,  MS 

NRL  Code  .58(X)  Wa.shington.  DC;  Code  5843  (F.  Rosenthal)  Washington.  DC;  Code  8441  (R.A,  Skop), 
Washington.  DC 


USCG  Code  2511  (Civil  Engrg).  Washington,  DC 

NSC  Cheatham  Annex.  PW’O.  W'illianishurg,  \’A;  Code  54  1,  Norfolk.  \'A:  Code  7(MI  Norfolk.  V'A:  SCE, 
Charleston.  SC:  SCE.  Norfolk,  \'A.  Seeurity  Offr  (Code  44).  Oakland.  CA 
NSD  SCE,  Subie  Bay.  RP 
CBC  4111,  OICC.  Great  Lakes.  IL 

NCC  LEAR  REGULATORY  COMMISSION  L  C  ,  Johnson,  Washington.  DC 

Nl'SC  DET  Code  .J,J22  (Varley)  New  London,  CT;  Code  EAI23  (R.S.  Munn).  New  London.  Cl;  Code  SB  .3,31 
(Brown).  Newport  RL  Code  lAlll  (G.  De  la  Cruz).  New  London  CT 
OFFICE  SECREI  ARY  OF  DEFENSE  ASD  IMRA&L)  Code  C'SSC'C  Washington.  DC;  OASD  (MRA&L) 

Dir.  of  Energy.  Pentagon,  Washington.  DC 
CNR  Code  4S1.  Bay  St.  Louis.  MS;  DEI'.  Dir,  Boston.  MA 
OCNR  Code  421  (Code  E  A.  Silva),  Arlington.  VA;  C'ode  7(K)F.  Arlington.  VA 
PERRY  OCEAN  ENG  R.  Pellen.  Riviera  Beach.  FL 

PHIBCB  1.  CO  San  Diego,  CA;  1.  ELCAS  Offer.  San  Diego,  Ca;  I.  P&E.  San  Diego.  CA;  2.  Co.  Norfolk.  VA 
PMTC  Code  5i).54-S.  Point  Mugu,  CA 

PWC  CO  (Code  61.3),  San  Diego,  CA;  Code  10,  Great  Lakes.  IL;  Code  10,  Oakland.  CA;  Code  101  (Library). 
Oakland.  CA;  Code  102,  Maim  Plan  &  Inspee,  Oakland,  CA;  Code  12.3-C.  San  Diego.  CA;  Code  2(K). 

Guam,  Mariana  Islands;  Code  .30V.  Norfolk,  VA:  Code  4(K(,  Pearl  Flarbor,  HI;  Code  400.  San  Diego.  CA; 
Code  420.  Great  Lakes.  IL;  Code  420.  Oakland.  CA;  Code  424.  Norfolk,  VA;  Code  425  (L.N.  Kaya,  P.E.). 
Pearl  Harbor,  HI:  Code  5()5A,  Oakland,  CA;  Code  590,  San  Diego.  CA;  Code  614,  San  Diego.  CA;  Code 
7(X).  San  Diego,  CA;  Dir  Maim  Dept  (Code  500),  Great  Lakes.  IL;  Dir.  Serv  Dept  (Code  4(X)).  Great 
Lakes.  IL;  Dir.  Transp  Dept  (Code  7(KI|,  Great  Lakes.  IL;  Dir.  Util  Dept  (Code  6(X)).  Great  Lakes.  IL; 
Library  (Code  1.34).  Pearl  Harbor,  HI;  Library.  Guam,  Mariana  Islands;  Library.  Norfolk,  VA;  Library, 
Pensaeola,  FL,  Library,  Yokosuka  JA.  Prod  Offr.  Norfolk.  VA;  Teeh  Library.  Subie  Bay,  RP;  Util  Offr, 
Guam,  Mariana  Island 
SEAL  TEAM  6,  Norfolk,  VA 
SPCC  PWO  (Code  08X).  Mechanicsburg.  PA 
SUPSHIP  Tech  Library.  Newport  News,  VA 
HAYNES  &  ASSOC  H.  Haynes.  P.E.,  Oakland.  CA 
UCT  ONE  OIC,  Norfolk.  VA 
UCT  TWO  OIC,  Port  Hueneme  CA 

U  S.  MERCHANT  MARINE  ACADEMY  Reprint  Custodian,  Kings  Point.  NY 

US  DEPT  OF  INTERIOR  Bur  of  Land  Mgmnt  Code  .58.3.  Washington  DC;  Natl  Park  Serv  (RMR/PC)  Denver. 
CO  80225 

US  GEOLOGICAL  SURVEY  Off.  Marine  Geology.  Piteleki.  Rcston  VA 

US  NATIONAL  MARINE  FISHERIES  SERVICE  Sandy  Hook  Lab.  Lib.  Highlands.  NY 

USCG  G-EOE-2/61  (Espinshade).  Washington,  DC;  G-EOE-4  (T  Dowd).  Washington.  DC;  Gulf  Strike  Team. 

Bay  St.  Louis.  MS:  LANT  Strike  Team.  Elizabeth  City,  NC;  Library  Hqirs,  Washington,  DC;  Pac  Strike 
Team,  Hamilton  AFB.  CA 

USCG  R&D  CENTER  CO,  Groton,  CT;  D.  Motherway.  Groton  CT;  Library,  Groton,  CT;  S  Rosenberg, 

Groton,  CT 

USCINC  PAC,  Code  J44,  Camp  HM  Smith.  HI 
USDA  Ext  Service  (T.  Maher)  Washington,  DC 

USNA  Ch.  Mech.  Engr.  Dept  Annapolis  MD;  Mech.  Engr.  Dept.  (Hasson).  Annapolis,  MD;  Mgr.  Engrg,  Civil 
Specs  Br.  Annapolis,  MD;  PWO,  Annafrolis,  MD 
USS  FULTON  WPNS  Rep.  Offr  (W-3)  New  York,  NY 
WATER  POWER  RESOURCES  SERVICE  (Smoak)  Denver,  CO 
ADVANCED  TECHNOLOGY  F.  Moss,  Op  Cen  Camarillo,  CA 
AMERICAN  CONCRETE  INSTITUTE  Detroit  MI  (Library) 

BERKELEY  PW  Engr  Div,  Harrison,  Berkeley.  CA 

CALIF.  DEPT  OF  FISH  &  GAME  Long  Beach  CA  (Marine  Tech  Info  Ctr) 

CALIF.  DEPT  OF  NAVIGATION  &  OCEAN  DEV.  Sacramento,  CA  (G.  Armstrong) 

CALIF.  MARITIME  ACADEMY  Library.  Vallejo,  CA 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY  Pasadena  CA  (Keck  Ref.  Rm) 

CALIFORNIA  STATE  UNIVERSITY  (Yen)  Long  Beach.  CA;  C.V.  Chelapati.  Long  Beach,  CA;  Dr.  Y.C. 

Kim.  Los  Angeles,  CA 

CITY  OF  AUSTIN  Resource  Mgmt  Dept  (G.  Arnold),Austin,  TX 
CITY  OF  LIVERMORE  Project  Engr  (Dackins),  Livermore.  Ca 
CLARKSON  COLL  OF  TECH  G.  Batson,  Potsdam  NY 

COLORADO  SCHOOL  OF  MINES  Dept  of  Engrg  (J.S.  Chung.  riiD)  Golden.  CO 

COLORADO  STATE  UNIVERSITY  Civil  Engr  Dept  (J.  Nelson)  f  Collins,  CO;  Civil  Engr  Dept  (W.A. 
Charlie)  Fort  Collins.  CO 

CORNELL  UNIVERSITY  Civil  &  Environ  Engrg  (F.  Kulhway),  Ithaca.  NY;  Library,  Ser  Dept,  Ithaca,  NY 

DAMES  &  MOORE  LIBRARY  Los  Angeles,  CA 

DUKE  UNIV  MEDICAL  CENTER  B.  Muga,  Durham  NC 

UNIVERSITY  OF  DELAWARE  (Dr.  S.  Dexter)  Lewes.  DE 

FLORIDA  ATLANTIC  UNIVERSITY  Boca  Raton  FL  (W,  Hartt);  Boca  Raton,  FL  (McAllister) 
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FLORIDA  lEC'HNOLOOlCAL  UNIVERSITY  Dr.  E.  Kahijian.  Mclhournc,  FL 
GEORC.IA  INSnn'TE  OF  TECHNOLOGY  Atlanta  C.A  (School  of  Civil  Engr..  Kahn) 

HARVARD  UNIVERSITY  Arch  Dept  (Mk  Kim).  Cambridge.  MA 
GEORC.IA  INSTHUTE  OF  TECHNOLOGY  Atlanta  GA  (B.  Mazanli) 

INSTFITriE  OF  MARINE  SCIENCES  Morchcad  City  NC  (Director) 

IOWA  STATE  LINIVERSIFY  Ames  lA  (CE  Dept.  Handy) 

WOODS  HOLE  OCEANOGRAPHIC  INST  Proj  Engr.  Woods  Hole.  MA 

LEHIGH  LINlVERSITY  Fritz  Engrg  Lab.  (Bcedic).  Bethlehem.  PA;  Linderman  Libr.  Ser  Cataloguer. 

Bcthicitem.  PA:  Marine  Geoleeh  Lab  (A.  Richards).  Bethlehem.  PA 
MAINE  MARHIME  ACADEMY  CASTINE.  ME  (LIBRARY) 

MICHIGAN  FECHNOLOGICAL  UNIVERSITY  Houghton.  Ml  (Haas) 

MIT  Engrg  Lib.  Cambridge.  MA;  Hydrodynamics  Lab  (Haricman).  Cambridge.  MA;  Library.  Cambridge.  MA; 

RV  Whitman.  Cambridge.  MA 
NATURAL  ENERGY  LAB  Library.  Honolulu.  HI 
NEW  MEXICO  SOLAR  ENERGY  INST.  Dr.  Zwibel  Las  Cruces  NM 
NYS  ENERGY  OFFICE  Library.  Albany  NY 

OREGON  SI  ATE  UNIVERSITY  (CE  Dept  Grace)  Corvallis.  OR;  CORVALLIS.  OR  (CE  DEPI .  BELL); 
C'orvalis  OR  (School  of  Oceanography) 

PENNSYLVANIA  STATE  UNIVERSITY  STATE  COLLEGE.  PA  (SNYDER);  State  College  PA  (Applied 
Rsch  Lab);  UNIVERSITY  PARK.  PA  (GOTOLSKI) 

PORT  SAN  DIEGO  Proj  Engr,  Port  Fac.  San  Diego,  CA 
PORTLAND  STATE  UNIVERSITY  H.  Migliore  Portland,  OR 

PURDUE  UNIVERSITY  Lafayette  IN  (Leonards);  Lafayette,  IN  (AltschaeffI);  Lafayette,  IN  (CE  Engr  Lib) 
SAN  DIEGO  .STATE  UNIV.  Dr  Krishnamoorthy,  San  Diego  CA;  I.  Noorany,  San  Diego,  CA 
SCRIPTS  INSTITUTE  OF  OCEANOGRAPHY  Deep  Sea  Drill  Proj  (Adams).  La  Jolla,  CA 
SEATTLE  U  Prof  Schwaegler  Seattle  WA 

SOUTHWEST  RSCH  INST  King.  San  Antonio.  TX;  R.  DeHart,  San  Antonio  TX;  San  Antonio.  TX 
STATE  UNIV.  OF  NEW  YORK  Buffalo,  NY;  Fort  Schuyler,  NY  (Longobardi) 

TEXAS  A&M  UNIVERSITY  College  Station  TX  (CE  Dept.  Herbich);  J.M.  Niedzwecki,  College  Station.  TX; 

W.B  Ledbetter  College  Station,  TX 
TEXAS  TECH  UNIVERSITY  Dept  of  IE  (Prof.  Ayoub),  Lubbock  TX 

UNIVERSITY  OF  ALASKA  Doc  Collections  Fairbanks,  AK;  Marine  Science  Inst.  College.  AK 
UNIVERSITY  OF  CALIFORNIA  A-(131  (Storms)  La  Jolla.  CA;  CE  Dept  (FX-IO).  Seattle,  WA;  CE  Dept 
(Taylor),  Davis,  CA;  La  Jolla  CA  (Acq.  Dept,  Lib.  C-()75A);  Naval  Arch  Dept,  Berkeley,  CA;  Prof  B.C. 
Gerwick,  Berkeley,  CA;  Prof  E.A.  Pearson,  Berkeley.  CA;  Prof  J.K.  Mitchell,  Berkeley.  CA 
UNIVERSITY  OF  CONNECTICUT  Library.  Groton,  CT 
UNIVERSITY  OF  DELAWARE  CivI  Engrg  Dept  (Chesson),  Newark,  DE 
UNIVERSITY  OF  FLORIDA  Florida  Sea  Grant  (C.  Jones),  Gainesville,  FL 
UNIVERSITY  OF  HAWAII  Library  (Sci  &  Tech  Div).  Honolulu,  HI 

UNIVERSITY  OF  ILLINOIS  Civil  Engrg  Dept  (Hall),  Urbana,  IL;  Library.  Urbana,  IL;  M.T.  Davisson. 
Urbana.  IL:  Metz  Ref  Rm,  Urbana,  IL 

UNIVERSITY  OF  MASSACHUSETTS  (Heronemus),  ME  Dept.  Amherst,  MA 
UNIVERSITY  OF  MICHIGAN  Ann  Arbor  MI  (Richart) 

UNIVERSITY  OF  NEBRASKA-LINCOLN  Lincoln,  NE  (Ross  Ice  Shelf  Proj.) 

UNIVERSITY  OF  NEW  HAMPSHIRE  P.  LaVoie,  Durham,  NH 
UNIVERSITY  OF  NOTRE  DAME  Katona,  Notre  Dame,  IN 

UNIVERSITY  OF  PENNSYLVANIA  Dept  of  Arch  (P.  McCleary).  Philadelphia.  PA;  Schl  of  Engrg  &  Applied 
Sci  (Roll),  Philadelphia,  PA 

UNIVERSITY  OF  RHODE  ISLAND  Narragansett  RI  (Pell  Marine  Set.  Lib.);  Wm,  D.  Kovacs.  Kingston,  R1 
UNIVERSITY  OF  SO.  CALIFORNIA  Hancock  Library,  Los  Angeles.  CA 
UNIVERSITY  OF  TEXAS  Inst.  Marine  Sci  (Library),  Port  Arkansas  TX 
UNIVERSITY  OF  TEXAS  AT  AUSTIN  (R.  Olson).  Dept  Civil  Engrg 

UNIVERSITY  OF  WASHINGTON  Applied  Physics  Lab,  Seattle.  WA.  CE  Dept.  (FX-10),  Seattle,  WA;  CE 
Dept,  Seattle,  WA;  Dept  of  Civil  Engr  (Dr.  Mattock).  Seattle  WA;  Library,  Seattle,  WA;  Pac  Marine 
Environ  Lab,  Seattle.  WA 

UNIVERSITY  OF  WISCONSIN  Great  Lakes  Studies.  Clr,  Milwaukee.  Wl 
VIRGINIA  INST.  OF  MARINE  SCI.  Library.  Gloucester  Point,  VA 
WESTERN  ARCHEOLOGICAL  CENTER  Library,  Tucson  AZ 
WOODS  HOLE  OCEANOGRAPHIC  INST  Doc  Lib.  Woods  Hole.  MA 
ALFRED  A.  YEE  &  ASSOC.  Librarian,  Honolulu.  HI 
AMETEK  Offshore  Res.  &  Engr  Div 
APPLIED  SYSTEMS  R.  Smith,  Agana.  Guam 
ARCAIR  CO.  D.  Young.  Lancaster  OH 
ARVID  GRANT  Olympia.  WA 

ATLANTIC  RICHFIELD  CO.  Engr  Serv  Grp  (J  Machemehi)  Dallas.  TX;  R  E.  Smith,  Dallas.  TX 
AUSTRALIA  Embassy  of  (Transport)  Washington.  DC;  Sch  Civil  Engr  &  Mining  (Prof  H.G.  Poulos),  Univ  of 
Sydney 


AWWA  RSC'H  FOUNDAI'ION  R.  Heaton,  Denver  CO 

BAITELLE-COLUMBUS  l.ABS  (D.  Frink)  Columbus,  OH;  (D.  Hackman)  Columbus,  OH 

BETHLEHEM  STEEL  CO.  Engrg  Dept  (Dismuke),  Bethlehem,  PA 

BRAND  INDUS  SERV  INC.  J.  Buehler,  Hacienda  Heights  CA 

BRITISH  EMBASSY  Sei  &  Tech  Dept  (Wilkins),  Washington.  DC 

BROWN  &  ROOT  Houston  TX  (D,  Ward) 

C  HEVRON  OIL  FIELD  RESEARCH  CO.  1.  Brooks.  La  Habra.  CA 
COLUMBIA  GULF  TRANSMISSION  CO.  Engrg  Lib.  Houston.  TX 
CONCRETE  TECHNOLOGY  CORP.  A.  Anderson,  Tacoma,  WA 
CONSTRUCTION  TECH  LAB  A.E.  Fiorato,  Skokie,  IL 
CROWLEY  MARITIME  SALVAGE  INC.  (B  Frost),  Willaimsburg.  VA 
DILLINGHAM  PRECAST  F.  McHale.  Honolulu  HI 
DIXIE  DIVING  CENTER  Decatur.  GA 
DRAVO  CORP  Pittsburgh  PA  (Wright) 

DURLACH,  O  NEAL,  JENKINS  &  ASSOC.  Columbia  SC 

EASTPORT  INTERNATIONAL  INC.  (J.H  Osborn)  Mgr,  West  Dtv,  Ventura.  CA 
ENERCOMP  H.  Amistadi,  Brunswick.  ME 
EVALUATION  ASSOC.  INC  MA  Fedele,  King  ol  Prussia.  PA 
EXXON  PRODUCTION  RESEARCH  CO  Houston.  TX  (Chao) 

FURGO  INC.  Library,  Houston,  TX 

GEOTECHNICAL  ENGINEERS  INC.  (R  F.  Murdock)  Principal.  W^inchestcr.  MA 

GOCJDYEAR  AEROSPACE  CORP  D/490,C2  (F  J  Stimler),  Akron.  OH 

GOULD  INC.  Tech  Lib,  Ches  Instru  Div  Glen  Burnie  MD 

GRUMMAN  AEROSPACE  CORP.  Tech  Info  Ctr.  Bethpage.  NY 

HALEY  &  ALDRICH,  INC.  HP  Aldrich.  Jr,  Cambridge.  MA 

KATSURA,  Y.  Coasult  Engr,  Ventura.  CA 

LAMONT-DOHERTY  GEOLOGICAL  OBSERVATORY  Palisades  NY  (McCoy) 

LIN  OFFSHORE  ENGRG  P.  Chow,  San  Francisco  CA 
LINDA  HALL  LIBRARY  Doc  Dept,  Kansas  City,  MO 
MARATHON  OIL  CO  Houston  TX 

MARINE  CONCRETE  STRUCTURES  INC.  W.A.  Ingraham.  Metairie.  LA 

MCDONNELL  AIRCRAFT  CO.  Sr  Engr.  Logistics,  St  Louis,  MO 

MOBIL  R  &  D  CORP  Offshore  Eng  Library,  Dallas,  TX 

MOFFATT  &  N/CHOL  ENGINEERS  (R.  Palmer)  Long  Beach.  CA 

EDWARD  K.  NODA  &  ASSOC  Honolulu,  HI 

PACIFIC  MARINE  TECHNOLOGY  (M.  Wagner)  Duvall,  WA 

PHELPS  ASSOC  P.A.  Phelps,  Rheem  Valley.  CA 

PORTLAND  CEMENT  ASSOC.  SKOKIE,  IL  (CORLEY;  SKOKIE.  IL  (KLIEGER);  Skokie  II,  (Rsch  &  Dev 
Lab.  Lib.) 

R  J  BROWN  ASSOC  (R.  Perera),  Houston,  TX 

RAYMOND  INTERNATIONAL  INC.  E  Code  Soil  Tech  Dept,  Pennsauken,  NJ;  J,  Welsh  Soiltech  Dept, 
Pennsauken.  NJ 

SANDIA  LABORATORIES  Library  Div.,  Livermore  CA;  Seabed  Progress  Div  4536  (D.  Talbert)  Albuquerque 
NM 

SCHUPACK  SUAREZ  ENGRS  INC  M.  Schupack,  South  Norwalk,  CT 
SEATECH  CORP.  MIAMI.  FL  (PERONI) 

SHANNON  <&  WILLSON  INC,  Librarian  Seattle.  WA 

SHELL  DEVELOPMENT  CO.  Houston  TX  (C.  Sellars  Jr.);  Houston  TX  (E.  Doyle) 

SHELL  OIL  CO.  E  &  P  CE,  Houston,  TX;  I.  Boaz,  Houston  TX 
SIMPSON  GUMPERTZ  &  HEGER  INC  Consulting  Engrs  (E.  Hill),  Arlington,  MA 
TANDEMLOC  INC  (J.  DiMartino  Jr)  Bayporl,  NY 
TIDEWATER  CONSTR.  CO  J  Fowler.  Virginia  Beach,  VA 

TRW  SYSTEMS  Engr  Library,  Cleveland.  OH;  M/S:  951/224  (P.K.  Dai).  San  Bernardino.  CA 
WESTINGHOUSE  ELECTRIC  CORP.  Annapolis  MD  (Oceanic  Div  Lib,  Bryan) 

WESTINSTRUCORP  Egerton,  Ventura,  CA 

WISS.  JANNEY,  ELSTNER.  &  ASSOC  Northbrook,  IL  (D.W.  Pfeifer) 

WM  CLAPP  LABS  -  BATTELLE  Library.  Duxbury,  MA 
WM  WOOD  &  ASSOC.  (D  Wood)  Metaire,  LA 

WOODWARD-CLYDE  CONSULTANTS  (R.  Cross).  Walnut  Creek,  CA;  Library.  West.  Reg.,  Walnut  Creek, 
CA 

ANTON  TEDESKO  Bronxville  NY 

BARA.  JOHN  P.  Lakewood.  CO 

BARTZ.  J  Santa  Barbara.  CA 

BRADFORD  ROOFING  T.  Ryan,  Billings,  MT 

BULLOCK  La  Canada 

F  HEUZE  Alamo.  CA 

F  W.  MC  COY  Woods  Hole.  MA 
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BEN  C.  GERWICK.  INC  San  Francisco,  CA 
HAYNES,  B.  Round  Rock,  TX 
LAYTON  Redmond,  WA 
CAPT  MURPHY  Sunnyvale.  CA 

MARINE  RESOURCES  DEV  FOUNDATION  N  T.  Monney.  Annapolis.  MD 

OSBORN.  JAS.  H.  Ventura,  CA 

PAULI  Silver  Spring,  MD 

PETERSEN,  CAPT  N.W.  Camarillo,  CA 

R.F.  BESIER  CE.  Old  Sayhrook.  CT 
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SUBJECT  CATEGORIES 

1  SHORE  FACIUTIES 

2  Construction  methods  ar>d  materials  (includirtg  corrosion 

control,  coatin9s) 

3  Waterfront  structures  (maintenance/deterioration  control) 

4  Utilities  (irKludirtg  power  cor>ditionin9) 

5  Explosives  safety 

6  Construction  equipment  arnj  machinery 

7  Fire  prevention  and  control 

8  Antenna  technology 

9  Structural  analysis  and  design  (including  numerical  and 

computer  techniques) 

10  Protective  construction  (including  hardened  shelters. 

shock  and  vibration  studies) 

1 1  Soil/rock  mechanics 

13  BEQ 

14  Airfields  ar>d  pavements 

15  ADVANCED  BASE  AND  AMPHIBIOUS  FACILITIES 

16  Base  facilities  (irtcluding  shelters,  power  generation,  water  supplies) 

17  Expedient  roads/airfields/bndges 

18  Amphibious  operations  (including  breakwaters,  wave  forces) 

19  Over-the-Beach  operations  (including  contameniation, 

materiel  transfer,  lighterage  and  cranes) 

20  POL  storage,  transfer  and  distribution 
24  POLAR  ENGINEERING 

24  Same  as  Advanced  Base  and  Amphibious  Facilities, 
except  limited  to  cold-region  environments 


28  ENERGV/POWER  GENERATION 

29  Thermal  conservation  (thernr^l  engineering  of  buildings.  HVAC 

systems,  energy  loss  measurement,  power  generation) 

30  Controls  and  electrical  conservation  (electrical  systems. 

energy  monitoring  and  control  systems) 

31  Fuel  flexibility  (liquid  fuels,  coal  utiliaation.  energy 

from  solid  waste! 

32  Alternate  energy  source  /geothermal  power,  photovoltaic 

power  systems,  solar  systems,  wir>d  systems,  energy  storage 
systems) 

33  Site  data  and  systems  integration  (energy  resource  data,  energy 

consumption  data,  integrating  energy  systems) 

34  ENVIRONMENTAL  PROTECTION 

35  Solid  waste  management 

36  Haaardous/toxic  materials  management 

37  Wastewater  management  and  sanitary  engineering 

38  Oil  pollution  removal  and  recovery 

39  Air  pollution 

40  Noise  abatement 

44  OCEAN  ENGINEERING 

45  Seafloor  soils  and  foundations 

46  Seafloor  construction  systems  ar>d  operations  (including 

diver  and  manipulator  tools) 

47  Undersea  structures  and  nriatenals 

48  Anchors  and  moorings 

49  Undersea  power  systems,  electromechanical  cables. 

and  connectors 

50  Pressure  vessel  facilities 

51  Physical  environment  (including  site  surveying) 

52  Ocean-based  concrete  structures 

53  Hyperbaric  chambers 

54  Undersea  cable  dynamics 
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